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Abstract 

Nonlinear  Quantitative  Feedback  Theory  (QFT)  is  applied  to  the  design  of 
an  flight  Control  system  (FCS)  to  control  maneuvers  of  an  agile  aircraft.  To  do 
this,  maneuvers  are  chosen  that  are  desirable  during  within  visual  range  combat. 
Plant  transfer  functions  are  developed  through  the  use  of  system  identification,  using 
the  input /output  time  histories  of  the  maneuvers.  During  the  system  identification 
process  the  qij{o)  are  found  directly  and  are  constrained  to  be  minimum  phase.  After 
identifying  the  plants,  the  cascaded  multiple-input  multiple-output  QFT  technique 
is  used  to  design  the  FCS.  During  the  design  of  the  FCS,  stability  boundaries  are 
generated  through  the  use  of  a  MatriXx  program.  To  prevent  the  formation  of  large 
RHP  poles,  F'boundries  are  developed.  Before  shaping  the  loops  it  is  found  that 
FCS  scheduling  is  needed  because  of  large  uncertainty  in  the  magnitude  and  phase 
of  the  plants.  Because  the  orders  of  the  designed  controllers  are  too  large  for  the 
YF-16  simulation  program,  reduced  order  controllers  are  developed  through  the  use 
of  straight  line  approximations  on  the  Bode  plot.  After  implementing  the  reduced 
order  controllers  in  the  nonlinear  simulation  program,  problems  developed  in  the 
modeling  process  are  encountered  that  prevent  the  simulation  of  the  FCS’s. 
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APPLICATION  OF  NONLINEAR  QFT 
TO  FLIGHT  CONTROL  DESIGN  FOR 
HIGH  ANGLE  OF  ATTACK  MANEUVERS 
WITH  THRUST  VECTORING 


I.  Introduction 


1.1  Background 

The  United  States  Air  Force  is  interested  in  weapon  systems  and  tactics  that 
will  give  its  pilots  an  advantage  over  the  enemy.  Because  of  relatively  new  weapons 
systems  such  as  the  all-aspect  missile,  tactics  can  be  developed  for  aircraft  that  will 
add  to  this  advantage.  If  aircraft  can  be  made  more  agile,  this  advantage  over  the 
enemy  can  be  improved.  Since  agile  aircraft  will  operate  in  highly  nonlinear  regions  of 
the  flight  envelope,  and  since  controlling  aircraft  in  these  regions  of  flight  will  depend 
on  very  fast  responses,  sophisticated  flight  control  systems  (FCS)  must  be  developed 
to  aid  the  pilot  in  controlling  these  aircraft.  This  thesis  is  part  of  a  continuing  effort 
to  improve  the  design  of  FCS  through  the  use  of  nonlinear  Quantitative  Feedback 
Theory  (QFT)  technique.  QFT  is  a  control  system  design  method  developed  by 
Dr.  Isaac  Horowitz  that  quantifies  the  amount  of  uncertainty  in  the  plant  and  uses 
this  information  to  reduce  the  negative  effects  of  using  feedback  in  the  controller. 
Nonlinear  QFT  uses  the  idea  of  quantifying  uncertainty  in  the  plant  and  extends  it 
to  quantifying  the  nonlinear  behavior.  By  quantifying  the  nonlinear  behavior  and 
the  uncertainty,  the  nonlinear  plant  can  be  rigorously  represented  in  this  design 
technique.  Nonlinear  QFT  was  used  in  previous  theses  by  Kobylarz  [12]  and  Miller 
[18].  Kobylarz  applied  nonlinear  QFT  with  pilot  compensation  to  the  design  of  a  FCS 
for  a  single-input  single-output  (SISO)  model  of  the  YF-16  aircraft.  Miller  applied 
the  design  techniques  to  the  design  of  a  FCS  for  a  multiple-input  multiple-output 
(MIMO)  model  of  the  YF-16,  where  there  were  two  inputs  and  two  outputs. This 
thesis  applies  nonlinear  QFT  to  design  a  FCS  for  a  three-input  three-output  model  of 
the  YF-16,  with  thrust  vectoring,  as  it  is  flown  through  highly  nonlinear  maneuvers. 
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1.2  Problem  Statement 

Since  agile  aircraft  will  be  able  to  maneuver  in  highly  nonlinear  regions  of  flight, 
a  method  of  flight  control  design  is  needed  to  give  an  acceptable  first-cut  design  of  a 
PCS  that  enables  the  pilot  to  control  the  aircraft  through  these  maneuvers.  In  the 
past  a  PCS  design  has  been  a  matter  of  cut  and  try.  By  using  a  method  such  as 
nonlinear  QPT,  the  amount  of  cut  and  try  can  be  reduced  significantly  [9]. 

1.3  Assumptions 

The  following  assumptions  are  made  for  this  thesis. 

•  Only  the  desired  outputs  are  of  interest  for  final  performance. 

•  The  robustness  of  MatriXx  transfer  function  algorithms  is  adequate  for  the 

design  process. 

The  first  one  is  the  basic  assumption  required  in  applying  nonlinear  QPT.  The 
second  is  made  because  Matrix*  has  limits  to  the  robustness  of  its  transfer  function 
algorithms. 

1.4  Scope 

This  thesis  applies  nonlinear  QPT  technique  to  the  design  of  a  PCS  for  the 
roll  (p),  pitch  (q),  and  yaw  (r)  channels  of  a  YP-16  simulation  program.  The  YP-16 
simulation  is  a  full  six-degrees-of-freedom  nonlinear  PORTRAN  simulation.  This 
simulation  includes  thrust  vectoring  in  the  pitch  and  yaw  channels.  The  PCS  is  de¬ 
signed  to  control  the  aircraft  through  two  specified  maneuvers.  The  maneuvers  are 
based  on  theoretical  maneuvers  for  agile  aircraft.  A  system  identification  method¬ 
ology  is  discussed  that  constrains  the  plants  to  be  minimum-phase  (MP).  QPT 
stability  boundaries  and  other  boundaries  are  computer  generated  to  aid  in  the  PCS 
design.  The  PCS  is  designed  to  have  the  time  responses  to  p,  q,  and  r  commands 
specified  in  MIL--STD-1707A  [17].  Order  reduction  on  the  controllers  is  accom¬ 
plished.  The  reduced  order  controllers  are  coded  in  PORTRAN,  and  installed  in 
the  YP-16  simulation  program.  Pinally,  problems  with  model,  that  become  evident 
during  the  nonlinear  simulation,  that  prevent  a  successful  simulation  are  discussed. 
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1.5  Standards 


The  thesis  sponsor  has  provided  specifications  on  the  response  to  control  inputs 
based  on  guidelines  presented  in  MIL-STD-1797A  [17].  The  specifications  on  the 
particular  maneuvers  were  developed  by  the  sponsor  [2]  [4]. 

1.6  Approach 

Maneuvers  are  developed  by  trial  and  error,  for  a  six  degree  of  freedom  non¬ 
linear  YF-16  simulation.  Input/output  time  histories  are  collected  that  adequately 
describe  the  maneuvers.  System  Identification  is  used  to  find  linear-time-invariant 
(LTI)  plants  based  on  these  time  histories.  Nonlinear  QFT  is  used  to  design  multi¬ 
ple  unity  feedback  FCSs  that  are  to  be  used  to  control  the  simulation  through  the 
prescribed  maneuvers.  Reduced  order  FCSs  are  coded  in  FORTRAN  and  scheduled 
based  on  the  maneuver.  The  actual  nonlinear  simulations  are  not  accomplished  due 
to  problems  that  arise  due  to  mismodeling. 

l.^  Summary 

This  thesis  consists  of  eight  chapters.  The  first  chapter  introduces  the  the¬ 
sis  problem.  The  second  chapter  gives  a  basic  overview  of  the  cascaded  nonlinear 
MIMO  QFT  method.  Chapter  3  explains  how  the  maneuvers  are  developed  and  also 
describes  the  YF-16  simulation.  Chapter  four  presents  the  manner  with  which  LTI 
plants  are  derived  from  input/output  time  histories.  The  fifth  chapter  describes  the 
method  used  in  obtaining  boundaries  on  the  forward  loop  transmission,  1.  Described 
in  Chapter  6  are  loop  shaping  and  the  design  of  the  prefilters.  The  seventh  chapter 
covers  the  problems  that  appear  during  the  nonlinear  simulations.  The  final  chapter, 
Chapter  8,  contains  a  summary,  conclusions,  and  recommendations. 
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11.  Cascaded  MEMO  Nonlinear  QFT 

2.1  Introduction 

The  essence  of  the  robustness  of  QFT  is  the  method’s  ability  to  quantify  uncer¬ 
tainty  [10].  This  chapter  outlines  the  QFT  iiiethod  and  how  nonlinear  plants  can  be 
modeled  as  LTI  plants  with  uncertainty.  First  the  QFT  method  is  reviewed.  Next, 
the  cascaded  MIMO  QFT  method  of  obtaining  MISO  effective  loops  is  explained. 
Finally,  the  manner  in  which  uncertain  LTI  plants  are  used  to  represent  nonlinear 
plants  is  presented. 

2.2  QFT  Method 

QFT  is  a  design  method  that  gives  the  designer  a  very  good  first  cut  at  the 
final  control  system  [9].  't  hrough  the  use  of  QFT,  one  can  design  a  fixed  parameter 
controller  based  on  a  range  of  operating  conditions.  Included  in  QFT  is  the  ability 
to  design  the  controller  to  cause  the  system  to  produce  outputs  that  meet  required 
specifications  based  on  particular  inputs.  The  overall  QFT  design  process  is  based 
on  the  diagram  presented  in  Figure  2.1,  where  the  elements  in  the  figure  can  be 
either  SISO  or  MIMO  [5]  [11].  Ucmd  represents  the  commanded  input  and  d\  ,  ^2 
are  disturbance  inputs.  G(s)  is  a  cascade  compensator  and  F(s)  is  a  prefilter.  In 
QFT  the  plant,  P(s)  ,  is  represented  by  a  set  of  J  LTI  transfer  functions,  represent¬ 
ing  the  plant  at  different  operating  conditions.  The  individual  transfer  functions, 
Pjis)  G  T’(s),  where  j  =  l,2,...,J,  are  chosen  by  the  designer  so  that  the  entire  opera¬ 
tional  envelope  of  the  plant  is  adequately  described.  The  QFT  technique  determines 
boundaries,  that  must  be  satisfied  by  a  nominal  loop  transmission  function,  which 
is  plotted  on  a  Nichols  chart  (NC).  The  boundaries  are  based  on  stability,  tracking, 
and  disturbance  rejection  performance  specifications.  In  this  thesis  only  stability 
bounds  are  specified  so  they  are  the  only  bounds  described  here.  The  nominal  loop 
transmission  (/^)  is  the  forward  loop  transmission  expressed  by: 

lo  =  GPo  (2.1) 

where  Po  represents  the  nominal  plant.  The  choice  of  Po  is  arbitrary,  but  a  particular 
Po  is  chosen  based  on  the  designer’s  preference  on  the  placement  and  shape  of  the 
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Figure  2.1.  QFT  Controller  Design 


boundaries.  To  form  the  stability  bounds,  templates  are  derived  based  on  V.  These 
templates  are  plots  of  the  magnitude  versus  phase,  for  a  given  frequency  a;;  for 
each  Pj  contained  within  V  =  {Pj}.  For  each  uji  a  stability  boundary  Bs{ju!)  can 
be  obtained  by  moving  the  template  around  the  M-contour  that  corresponds  to  the 
required  phase  margin  angle  7  ,  and  plotting  ^he  locus  of  points  described  by  P^,. 
Once  boundaries  are  obtained  for  a  set  of  specified  values  of  w,,  C  can  be  formed 
on  the  NC,  with  the  restraint  that  at  each  Wj,  COw,)  must  lie  on  or  above  Bs{juj). 
After  lo  is  formed,  G(s)  can  be  derived  from  Equation  (2.1).  Once  the  feedback  loop 
is  formed,  F(s)  can  be  designed.  A  modified  method  is  used  to  design  F(s),  This 
method  involves  the  following  steps: 

1.  A  contro'  ratio  is  obtained  for  each  Lj{s)  ~  G{s)Pj{s). 

2.  The  frequency  response  is  obtained  for  each  Tl^{s) 

3.  The  Bode  plots  of  all  and  the  frequency  responses  of  the  tracking  bounds, 
Tfi^  and  ,  are  ploUed  on  the  same  bode  plot. 

4.  F{s)j  is  designed  so  that  all  Lm[F(s)j7£,j]  lie  within  the  bounds  given  by 
h’FTvj]  and  Lm[Ti,J 

2.3  Cascaded  MIMO  QFT 

In  general,  the  MIMO  QFT  design  method  involves  transforming  the  system 
of  Figure  2.1,  with  no  external  disturbances,  into  equivalent  multiple-input  single¬ 
output  (MISO)  loops  as  shown  in  Figure  2.2  for  the  3X3  case.  This  transformation 
is  justified  through  the  use  of  the  fixed  point  theorem  [11].  These  MISO  loops  are 
represented  by  equations  of  the  form  of  Equation  (2.2),  where  is  the  control  ratio 
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relating  the  output  to  the  input,  /,  is  the  prefilter,  /,  is  the  forward  loop 
transmission,  and  d,j  represents  a  disturbance  input. 


fijh  ~h 

^  +  h 


(2.2) 


The  q,j  in  these  equations  are  the  reciprocal  of  the  terms  contained  in  the  inverse 
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Figure  2.2.  3X3  MISO  Equivalent  Loops 


of  the  square  plant  matrix  P=  [pi^].  That  is  P*=[p*^j=F“*  where  \p*-]  =  and  Q 
=  The  disturbances  in  these  MISO  loops  represent  the  cross-coupling  between 
loops.  To  simplify  the  controller,  the  G(s)  and  F(s)  matrices  used  in  this  thesis  are 
diagonal,  but  non-diagonal  Ci(s)  and  F(s)  can  be  used.  The  design  of  the  first  loop 
is  accomplished  using  the  method  described  for  the  MISO  case.  After  the  first  loop 
is  designed,  the  set  of  effective  plants,  that  accounts  for  the  effect  of  the  first 

loop  on  the  succeeding  loops  is  derived  using  Equ  ations  (2.3)  through  (2.6)  [8]. 


.  _  fijh(k)  ~h 

1  +  kk) 


(2.3) 


h{k)  —  9iQij{k) 


(2.4) 
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_  9q(fc-i)(l  +  /t(fc-i)) 

‘•'(''I  1  -  -1-  /.(jt-i) 

(2.5) 

9»i(fc— l)9fcA:(fc— 1) 

(2.6) 

Where  k  in  these  equations  indicates  the  kth  loop  that  is  previously  designed.  These 
are  the  same  effective  plants  described  by  Houpis,  [11]  and  Horowitz,  [7],  but  these 
equations  lend  themselves  to  a  recursive  approach  to  finding  effective  plants.  This 
approach  is  used  here  because  the  number  of  functions  used  in  the  computer  algo¬ 
rithm  is  reduced  allowing  a  simpler  implementation.  For  the  second  loop  only  922(1) 
needs  to  be  calculated,  but  for  the  third  loop  933(1),  922(1)5  932(1)?  and  923(1)  all  need 
to  be  calculated.  With  922(1)?  the  second  loop  can  be  designed  in  a  similar  manner 
as  the  first  loop.  Once  the  second  loop  is  determined,  the  set  of  effective  plants  for 
the  third  loop  ^33(2)  is  calculated  based  on  Equations  (2.3)  through  (2.6)  with  j=3, 
i=3,  and  k=2.  Finally  the  third  loop  is  designed  in  the  manner  of  the  previous  loops. 
After  the  three  loops  are  determined,  the  diagonal  prefilters  are  designed. 


2.4  Nonlinear  QFT 

The  nonlinear  behavior  of  a  plant  is  quantified  using  the  same  method  that 
quantifies  the  uncertainty  of  a  plant.  That  is,  a  set  of  LTI  transfer  functions  {Vi), 
is  found  that  represents  a  nonlinear  plant  IF,  over  a  set  of  acceptable  outputs  O 
where  VFj  6  W  ,  Figure  2.3.  A  new  Vi  must  be  found  for  every  Wj  6  W.  It  hcis  been 
proven,  using  the  fixed  point  theorem,  that  an  acceptable  design  based  on  Vi  will 
satisfy  the  design  requirements  of  VW,  C?,  [9].  The  design  process  involves  first  finding 
V,  where  Vi  €  V.  To  find  a  particular  Vi  an  envelope  of  desired  responses  must  be 
specified  for  the  particular  operating  condition.  These  are  all  the  possible  responses 
that  are  required  from  the  plant  at  that  operating  condition.  A  set  of  outputs  (9, 
are  chosen  to  adequately  describe  Wi,  within  this  envelope.  The  number  of  outputs 
needed  will  depend  on  the  particular  W^.  The  set  of  inputs  Xi  that  are  required  to 
produce  C?,,  are  used  to  find  the  set  of  LTI  transfer  functions,  Vi.  Finding  Vi  can 
be  accomplished  by  either  analytical  means  or  as  in  the  case  of  this  thesis  system 
identification.  Next  a  set  of  templates  is  generated,  based  on  V,  as  in  the  SISO  case. 
From  this  point  the  design  method  is  the  same  as  the  LTI  design  method. 
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Figure  2.3.  Nonlinear  Diagram 


2.5  Conclusion 

This  chapter  reviews  the  basic  techniques  of  the  QFT  method.  A  basic  outline 
of  MIMO  QFT  is  given  and  the  cascaded  method  of  effective  plant  computation 
is  explained.  Finally,  the  difference  between  nonlinear  QFT  and  linear  QFT  is  ex¬ 
plained.  In  the  next  chapter  the  desired  maneuvers  and  the  agile  aircraft  used  in 
this  thesis  are  described. 
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III.  Desired  Maneuvers  and  Agile  Aircraft 


3.1  Introduction 

‘Simply  stated,  fighter  agility  is  that  characteristic  of  a  fighter  aircraft  which 
enables  it  to  outmaneuver  its  adversary  for  a  favorable  position  from  which  to  launch 
its  weapons  [21].’  Agility  is  an  important  part  of  air-to-air  combat  that  occurs  within 
visual  range.  Combat  will  be  conducted  within  visual  range  for  many  reasons  such  as 
large  enemy  numbers,  stealth  aircraft,  faulty  missiles,  and  visual  identification  [21] 
[6].  Since  an  agile  fighter  is  desirable  considerable  interest  exists  in  techniques  that 
will  improve  agility.  As  Nguyen  and  Gilbert  said,  ‘high  agility  requires  the  ability 
to  generate  substantial  controllable  angular  rates  over  an  expanded  angle  of  attack 
(AOA)  range  [i9].’  An  aircraft  capable  of  maneuvering  at  high  AOA  can  change 
its  attitude  very  quickly  and  therefore  point  its  nose  at  a  target  very  quickly  [21]. 
This  chapter  describes  the  development  of  the  maneuvers  that  the  PCS  is  designed 
to  control  and  the  aircraft  model  used  in  this  thesis. 

3.2  Maneuvers 

Since  increasing  the  AOA  of  an  aircraft  increases  drag,  maneuvers  at  high 
AOA  cause  a  large  loss  of  speed.  This  loss  of  speed  translates  into  a  loss  of  energy 
for  combat  maneuvering.  Loss  of  energy  can  put  an  attacking  aircraft  at  a  distinct 
disadvantage  if  the  first  shot  misses  or  if  there  are  multiple  targets  [21].  Tamrat 
suggests  that  the  most  important  use  of  high  AOA  maneuvers  is  in  a  defensive  roll. 
One  of  the  maneuvers  is  to  pitch  up  to  a  certain  angle  and  then  point  the  nose  at 
the  target  using  the  aircraft’s  yaw  axis.  Another  maneuver  is  to  increase  the  AOA 
quickly,  which  causes  a  rapid  decrease  in  the  aircraft’s  forward  airspeed,  causing  the 
pursuing  aircraft  to  move  in  front  of  the  defending  aircraft  [21].  To  regain  energy 
from  these  high  AOA  maneuvers,  the  method  is  to  pull  the  aircraft’s  nose  down 
quickly  to  unload  the  wings  and  rapidly  increcise  airspeed  [19].  Another  maneuver 
that  can  be  useful  is  a  banked  turn  at  a  high  AOA.  This  maneuver  would  give  the 
aircraft  the  ability  to  turn  toward  the  enemy  quickly  allowing  a  first  shot  [2]  [21]. 

3.3  Thrust  Vectoring 

Thrust  vectoring  is  an  attractive  control  to  use  on  a  fighter  aircraft  because 
it  ‘provides  control  moments  that  are  essentially  uncoupled  from  airframe  aerody- 
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namics  [3].’  At  flight  conditions  where  aerodynamic  control  surfaces  have  little  or 
no  effectiveness,  thrust  vectoring  provides  control  power.  Examples  of  these  flight 
conditions  are  low  airspeed  and  high  AO  A  [3].  Nguyen  and  Gilbert  state  that  thrust 
vectoring  matches  or  exceeds  the  control  effectiveness  needed  to  control  aircraft 
during  low  speed  flight  [19].  Aircraft  designers  must  use  larger  control  surfaces  to 
compensate  for  low  speed  flight.  Thrust  vectoring  allows  the  aircraft  designer  to 
make  these  surfaces  smaller,  thus  reducing  drag  and  weight  [3].  At  higher  AOA’s 
the  control  surfaces  don’t  provide  enough  control  authority  to  maneuver  the  air¬ 
craft.  During  these  flight  conditions  control  surfaces  can  be  augmented  or  replaced 
by  thrust  vectoring.  At  AOA’s  above  the  stall  AOA,  thrust  vectoring  must  be  used 
to  control  the  aircraft  [21].  During  their  research  on  three  different  fighter  configu¬ 
rations,  Capone  and  Mason  [3]  found  that  the  pitching  moment  coefficient  and  the 
lateral  aerodynamic  coefficients  varied  linearly  with  the  amount  of  thrust  vectoring. 
No  cross-coupling  of  the  thrust  vectoring  controls  existed  between  the  lateral  and 
longitudinal  forces  and  moments.  Thrust  vector  control  power,  for  their  design,  was 
a  function  of  Mach  number  and  nozzle  pressure  ratio  .  In  his  design  for  a  thrust  vec¬ 
toring  controller  of  the  F-15,  Anderson  [1]  connected  the  pitch  thrust  control  input 
to  the  same  command  as  the  collective  stabilators.  The  yaw  thrust  vectoring  was 
controlled  by  use  of  a  command  proportional  to  the  sideslip  angle.  Because  the  F-15 
has  two  engines  a  roll  thrust  vectoring  control  was  used  and  this  was  commanded  by 
the  differential  stabilator  command.  His  design  was  tested  in  a  flight  simulator  with 
USAF  F-15  pilots  and  with  simulated  pilot  commands.  Anderson’s  results  indicate 
that  his  thrust  vector  augmented  F-15  was  superior  for  most  tasks  when  compared 
to  the  unaugmented  F-15  model.  One  major  drawback  cited  was  the  loss  of  airspeed 
in  extended  maximum  turns. 

3.4  Aircraft  Model 

The  aircraft  model  used  in  this  thesis  is  a  FORTRAN  implemented  simulator 
of  the  YF-16.  This  is  a  nonlinear  simulation  and  represents  the  full  six  degrees 
of  freedom  equations  of  motion.  The  simulator  is  provided  by  Tom  Cord  of  FDL 
[4].  The  YF-16  is  a  prototype  version  of  the  F-16,  but  it  has  less  internal  avionics 
and  therefore  it  is  much  lighter  and  more  agile  [2].  Included  in  this  simulation  is  a 
mechanism  to  add  thrust  vectoring  to  the  YF-16.  Thrust  vectoring  is  modeled  as  a 
controllable  moment  driven  through  an  actuator,  that  is  produced  by  multiplying  a 
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portion  of  the  forward  thrust  by  the  distance  from  the  center  of  gravity  of  the  aircraft 
to  the  engine  nozzle,  in  the  pitch  and  yaw  directions,  see  Equations  (3.1)  through 
(3.3)  and  Figure  3.1,  where  STYpUch  and  STVyaw  are  thrust  deflection  angles,  is 
the  thrust  along  the  x  body  axis,  and  My  ,  M*  are  moments  produced  by  the  thrust 
deflection  about  the  y  and  z  body  axes. 

T 

Tx  =  -7===========  (3.1) 

^1  +  tan^{STVpitch)  +  tan'^  {STVyaw) 


My  =  dcgTxtan{6TVpitch) 
Mx  =  dcgTxtan{STVyaw) 


(3.2) 

(3.3) 


The  effect  of  decreased  forward  thrust  is  also  accounted  for  in  Equation  (3.1).  Cord’s 
simulation  program  contains  a  FCS  which  is  denoted  in  this  thesis  as  CFCS.  The 
CFCS  is  modified  to  a  version  described  in  Chapter  4,  shown  in  Figure  3.2,  and 
denoted  MFFCS  .  This  modification  is  done  in  order  to  obtain  the  LTI  models  for 
the  nonlinear  plant  as  described  in  Chapter  4. 

3.5  Maneuvers 

The  maneuvers  are  chosen  to  represent  the  types  of  maneuvers  that  Tamrat 
suggested.  Three  maneuvers  are  chosen  but,  due  to  time  constraints  only  two  ma¬ 
neuvers  are  modeled.  The  first  maneuver  is  a  quick  pitch-up  to  a  certain  pitch  angle, 
hold  that  angle  and  then  pitch-down  to  unload  the  wings  and  regain  energy.  The 
second  maneuver  consists  of  a  quick  pitch- up  and  then  a  yaw/roll  around  the  velocity 
vector.  The  third  maneuver  was  to  be  a  high  alpha  angle  turn.  To  simulate  the  ma¬ 
neuvers  the  original  control  system  is  modified  as  shown  in  Figure  3.2.  The  pitch  and 
roll  channel  FCSs  are  removed  and  replaced  by  simple  feedbacks.  The  yaw  channel 
FCS  is  left  intact.  These  modifications  are  made  to  remove  limiting  devices  in  the 
original  FCSs  that  restricted  the  AOA  and  angular  rates  of  the  aircraft  simulation. 
There  is  further  discussion  of  the  control  configuration  in  Chapter  4.  Inputs  to  the 
modified  FCS  are  in  the  form  of  vectors  which  are  indexed  in  0.05  sec  increments. 
These  vectors  are  developed  using  a  MatriXx  program  that  facilitates  trial  and  error 
input  development.  The  simulation  program  is  modified  to  read  in  the  input  vectors 
and  write  out  Matrix*  programs  that  cause  Matrix*  to  plot  the  response  verses  time 
of  several  variables  of  interest.  Before  developing  the  specifications  on  the  maneu¬ 
vers,  simulation  runs  are  made  to  learn  how  to  maneuver  the  simulated  aircraft  in  a 
manner  similar  to  the  desired  maneuvers.  This  is  necessary  to  investigate  how  the 
aircraft  reacts  to  different  sets  of  inputs  and  controller  gains.  After  each  simulation 
run,  angular  rates,  inputs  to  the  plant,  Euler  angles,  altitude,  speed,  AOA,  sideslip 
angle,  and  control  deflections  are  plotted  versus  the  time  of  the  simulation  as  shown 
in  the  example  in  Figures  (3.3)  through  (3.5)  and  Appendix  A.  From  these  plots 
one  can  develop  a  good  idea  of  what  happens  to  the  aircraft  during  the  maneuvers. 
At  this  point,  inputs  are  developed,  by  trial  and  error,  to  manipulate  the  simulated 
aircraft  in  order  to  fine  tune  the  gains  of  the  simplified  FCS.  One  problem  that  arises 
is  an  excess  /3  develops  during  some  of  the  simulations.  Because  of  stresses  developed 
on  the  airframe  due  to  /?,  this  angle  is  to  be  kept  below  three  degrees.  This  problem 
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is  corrected  by  adjusting  the  thrust  vector  gains  and  increasing  the  engine  thrust  [2] 
[4].  The  mechanisim  used  to  increase  the  thrust  from  trim  thrust  level  to  maximum 
thrust,  causes  a  step  increase  of  the  thrust  at  the  start  time  of  the  simulation.  This 
step  increase  in  thrust  is  a  control  input  to  the  plant  that,  because  of  an  oversight, 
is  not  modeled  in  the  development  of  the  FCS  and  causes  complications  in  latter 
stages  of  the  design. 

3.6  Specifications 

Since  there  are  no  guidelines  on  the  specifications  on  these  types  of  maneuvers, 
specifications  are  developed  that  are  based  on  maximum  and  minimum  attainable 
rates  and  accelerations  for  the  aircraft  at  certain  attitudes  and  the  experience  of 
Barfield  and  Cord  [2]  [4].  For  the  pitch-up  angular  acceleration  the  specification 
is:  the  angular  acceleration  available  from  eighty  percent  of  the  difference  between 
the  maximum  control  power  and  the  control  power  needed  for  stabilization.  The 
pitch-down  specification  is  given  by  the  angular  acceleration  based  on  the  stability 
derivative  Cmf  =0.05.  These  two  specifications  on  pitch  rate  acceleration  q  are 

*min 

calculated  with  Equations  (3.4)  through  (3.9). 
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(3.5) 

(3.6) 

(3.7) 

(3.8) 
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The  actual  values  for  the  variables  are  printed  out  from  the  simulation  program.  An 
additional  specification  is  made  to  keep  stresses  on  the  airframe  within  acceptable 
limits,  thus  the  sideslip  angle  ^  can  not  be  more  than  3  deg.  The  roll  portion  of 
the  maneuvers  is  to  accomplish  a  stability  axis  roll  at  three  different  rates.  Equation 
(3.10)  is  used  to  convert  the  body  axis  rates  to  stability  axis  rates.  The  specifica- 
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tions  for  the  maneuvers  are  contained  in  Table  3.1.  Consolidated  plots  of  all  of  the 
controlled  variables  and  the  velocity  vector  roll  are  shown  in  Figures  (3.6)  through 


(3.9). 


Pstab  =  pcos(a)  + 


ptan(a) + 


g  cos(0) 


sin(Q!) 


(3.10) 


Table  3.1.  Specifications  on  Maneuvers 


Variable 

Minimum 

Maximum 

Positive  q 

15  deg/sec 

50  deg/sec 

Negative  q 

-15  deg/sec 

-50  deg/sec 

<7 

15  deg/sec/sec 

50  deg/sec/sec 

Velocity  Vector  Roll 

30  deg/sec 

60  deg/sec 
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Figure  3.3.  Example  Maneuver  #211 
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Figure  3.5.  Example  Maneuver  #211,  cont’d 
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Figure  3.6.  Consolidated  Pitch-up/Pitch-down  Maneuvers 
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Figure  3.7.  Consolidated  Pitch-up/Pitch-down  Maneuvers,  cont’d 
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Figure  3.8.  Consolidated  Pitch-up/Pitch-down  Maneuvers,  cont’d 
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Figure  3.9.  Consolidated  Pitch-up/Roll  Maneuvers 
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Figure  3.11.  Consolidated  Pitch-up/Roll  Maneuvers,  cont’d 

3-16 


3. 7  Summary 

This  chapter  presents  the  background  for  the  choice  of  maneuvers  and  the 
aircraft.  The  maneuvers  are  chos?n  to  emulate  high  AO  A  maneuvers  that  are  theo¬ 
rized  advantages  for  future  aerial  combat.  The  use  of  thrust  vectoring  with  the  YF-16 
simulation  produces  a  fighter  aircraft  model  that  can  operate  in  an  expanded  flight 
envelope.  This  chapter  also  presents  the  specifications  on  the  maneuvers  and  hence 
the  FCS.  Chapter  4  describes  the  model  development  based  on  system  identification. 
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IV.  Model  Development 


4.1  Introduction 

This  chapter  describes  the  methods  used  and  discusses  some  of  the  obstacles 
encountered  during  development  of  the  mathematical  models  that  are  used  to  de¬ 
scribe  the  aircraft  in  the  nonlinear  QFT  design  process.  According  to  Ljung  and 
Soderstrom  developing  a  mathematical  model  of  a  physical  system  based  on  in¬ 
put/output  time  histories  is  called  system  identification  [14].  There  are  two  steps  to 
this  task:  first  is  deciding  on  a  family  of  candidate  models  and  second  is  selecting 
the  particular  model  from  that  set  that  adequately  relates  the  input /output  time 
histories.  The  adequacy  of  this  model  is  an  issue  the  developer  must  address  based 
on  the  application  intended  for  the  model.  There  are  two  major  classes  of  identifica¬ 
tion  techniques,  recursive  and  batch.  Recursive  techniques  are  predominantly  used 
in  on-line  applications  such  as  adaptive  controllers.  These  techniques  usually  oper¬ 
ate  on  one  segment  of  the  discrete  time  history  per  iteration.  Batch  identification 
techniques  operate  on  the  entire  time  history  at  one  time  and  are  used  for  identifying 
systems  where  the  entire  time  histories  are  available.  Because  all  of  the  time  history 
is  known,  a  batch  identification  technique  finds  an  optimized  model  of  the  system, 
compared  to  models  produced  by  recursive  techniques.  Since  all  of  the  time  history 
is  available  for  this  thesis,  a  batch  method  is  used. 

4.2  Previous  Research 

Kobylarz  and  Miller,  in  their  theses,  used  a  batch  technique  for  developing 
their  mathematical  models.  This  technique  was  based  on  a  technique  devised  by 
Golubev  [12]  [18].  Given  a  model  for  the  transfer  function  or  functions,  this  technique 
calculates  the  parameters  of  the  model  from  the  defining  mathematical  equations 
using  a  least  squares  approach.  A  good  explanation  of  this  technique  can  be  found 
in  Miller’s  thesis  [18].  Kobylarz  [12]  used  this  technique  on  a  SISO  problem  and 
Miller  extended  the  technique  for  his  MIMO  design. 

4.3  Other  Techniques 

Part  of  the  goal  of  this  thesis  is  to  explore  other  methods  of  system  identifi¬ 
cation.  The  majority  of  techniques  used  in  batch  identification  are  based  on  least 
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squares,  extended  Kalman  Filter,  and  maximum  likelihood.  Techniques  represent¬ 
ing  all  three  of  these  mathematical  theories  are  implemented  in  Matrixx  [16].  The 
least  squares  approach  uses  linear  algebra  to  manipulate  the  time  histories  and  then 
uses  least  squares  to  approximate  the  transfer  function  model  parameters.  By  using 
an  extended  Kalman  filter,  the  EKF  function  in  Matrixx  can  estimate  the  states  of 
the  syst'^m.  The  state  estimation  can  be  used  as  is  or  transformed  into  a  transfer 
function  form.  This  thesis  uses  a  MatriXx  batch  technique  based  on  a  maximum 
likelihood  function  called  maxlike.  There  axe  two  criteria  used  in  the  selection  of  the 
maxlike  function.  The  identification  technique  has  to  be  able  to  handle  MIMO  plant 
identification.  Because  of  numerical  accuracy  limitations  in  Matrixx  ,  the  function 
must  be  able  to  handle  continuous  transfer  function  matrices  in  the  Laplace  domain. 
The  maxlike  function  meets  both  criteria. 

4.4  Maximum  Likelihood  Function 

Maxlike  is  a  system  identification  function  based  on  a  maximum  likelihood 
function.  Essentially  what  maxlike  does  is  iteratively  execute  a  model,  defined  by 
the  user,  that  relates  the  output  time  history  data  to  the  input  time  history  data. 
This  model  can  be  any  MatriXx  executable  program.  During  the  iterations  maxlike 
uses  a  maximum  likelihood  function  to  estimate  parameters  contained  in  the  model 
to  obtain  a  best  fit  of  the  model  to  the  input/output  time  histories.  An  initial 
parameter  vector  is  supplied  by  the  user.  The  maxlike  function  iterates  until  it 
converges  on  a  solution  or  until  the  user  defined  maximum  number  of  iterations  are 
exhausted.  The  maxlike  function  executes  the  model  once  for  each  parameter  for 
each  point  in  time  so  that  the  number  of  parameters,  the  complexity  of  the  model 
and  the  number  of  time  points  all  have  an  effect  on  the  length  of  time  required  to 
find  a  solution. 

4.5  Testing  Maxlike 

To  become  familiar  with  and  test  the  maxlike  function  the  F-16  simulation 
program  Genesis  is  used  to  obtain  linearized  equations  of  motion  for  an  F-16  traveling 
at  0.9  Mach  at  30,000  ft.  Using  MatriXx  ,  step  responses  are  produced  for  the 
equations  and  then  using  a  state-space  model,  a  system  identification  run  is  made.  A 
state-space  model  is  used  because  the  linearized  equations  produced  by  Genesis  are  in 
state-space  form.  A  step  response  for  the  resulting  state-space  system  is  calculated 
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and  compared  to  the  original  step  response.  The  two  responses  fit  together  very 
closely  and  therefore  the  new  state  space  is  a  good  candidate  model  of  the  first. 
There  is  a  concern  with  cross-terms  in  the  resulting  equations  being  similar  to  those 
in  the  original  equations,  but  according  to  Horowitz,  in  using  the  nonlineaj  QFT 
method,  the  equivalent  system  must  reproduce  the  original  output  given  the  original 
input  and  no  constraints  are  made  as  to  the  appearance  of  the  equivalent  plant  [9]. 
During  this  testing,  it  is  apparent  that  the  quality  of  resulting  equivalent  systems 
is  a  function  of  the  initial  parameters,  the  form  of  the  model,  and  the  length  of 
time  of  the  data.  Furthermore  the  length  of  time  the  identification  process  takes 
is  a  function  of  the  computer  system  used.  For  one  particular  set  of  data  and  one 
particular  model  the  identification  process  takes  as  long  as  6  hours  on  a  VAX  to  30 
minutes  on  a  Sparc  workstation.  This  time  is  a  significant  factor  when  ones  considers 
that  a  laxge  number  of  equivalent  plants,  26  in  this  case,  are  necessary  to  model  a 
relatively  large  flight  envelope. 

4-6  Simulator  Compensator  Configuration 

After  verifying  that  maxlike  produces  acceptable  results  for  the  test  case,  the 
next  step  is  to  set  up  the  simulation  program  to  be  able  to  generate  the  necessary 
time  histories.  Cord’s  YF-16  simulator  program  is  modified  by  removing  the  original 
compensator  for  the  pitch  and  roll  channels  and  modifying  the  compensator  for  the 
yaw  channel.  As  shown  in  Figure  3.2,  the  original  pitch  and  yaw  compensators  are 
replaced  by  simple  feedback  loops  that  consist  of  negative  feedback  of  the  variable  of 
interest  and  the  insertion  of  a  gain  in  the  forward  loop.  This  type  of  compensation  is 
used  so  that  the  commanded  output  of  the  aircraft  is  a  signal  similar  to  the  desired 
output.  Another  reason  to  modify  the  existing  pitch  and  roll  compensation  is  the 
original  compensation  contained  devices  to  limit  the  angle  of  attack. 

4-7  Yaw  Channel  Configuration 

Because  the  goal  of  a  velocity  vector  roll  is  to  keep  the  sideslip  ^  at  a  minimum, 
some  kind  of  feedback  is  required.  This  feedback  can  be  made  part  of  the  plant  as 
shown  in  Figure  3.2.  In  this  figure  the  original  yaw  channel  compensation  is  used 
and  the  feedback  signals  representing  the  yaw  and  ^  are  separated.  /3  feedback  can 
also  be  accomplished  after  the  inner  loop  compensator  is  designed.  In  this  situation 
the  original  ^  feedback  is  used  to  obtain  the  desired  output.  Figure  4.1  shows  how  ^ 
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feedback  is  used  for  plant  generation  but  is  not  part  of  the  plant.  As  explained  later 
in  this  chapter,  there  is  a  big  difference  in  how  each  of  these  plants  can  be  modeled 
with  the  maxlike  function. 


actuators 

I - 1 

I  20  I 

f  s+20  1 


YF-16 


0  Feedback! 


Figure  4.1.  Simulation  Configuration  without  0  Feedback  included  in  the  Plant 


4-S  Thrust  Vectoring 

In  his  simulation  program,  Cord  designed  a  method  to  actuate  a  change  in  the 
direction  of  the  thrust  nozzle.  For  this  simulation,  based  on  Cords’s  experience  [4], 
actuators  similar  to  those  of  the  control  surfaces  are  used  for  the  thrust  vectoring. 
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The  control  signals  to  the  pitch  actuators  and  the  yaw  actuator  are  connected  to  the 
pitch  and  yaw  thrust  vectoring  actuators  through  gain  blocks  as  shown  in  Figure  4.1. 
By  connecting  the  thrust  vectoring  this  way  a  square  plant  matrix  is  achieved,  thus, 
the  task  of  developing  a  weighting  matrix  is  eliminated.  The  gains  in  the  thrust 
vectoring  loops  take  the  place  of  the  weighting  matrix  and  become  part  of  the  plant. 

4.9  Simulations 

After  the  simulation  program  is  modified  the  input/output  time  histories  are 
produced  based  on  the  specified  maneuver  envelopes.  To  produce  the  desired  output 
from  the  simulation  program,  input  command  time  history  vectors  are  produced 
using  MatriXx  .  To  find  the  desired  response  of  the  aircraft,  the  response  to  these 
input  vectors,  the  variables  of  interest  are  plotted  versus  time.  At  first,  simple  step 
commands  are  applied  to  the  inputs  and  the  gains  are  adjusted  to  obtain  the  best 
response.  The  best  response  is  one  that  follows  the  command  as  closely  a.s  possible 
without  causing  oscillations  in  the  control  surfaces.  After  many  trials  unity  gain  for 
pitch,  roll,  and  yaw  seem  to  give  best  results  for  all  maneuvers.  The  thrust  vectoring 
gains  have  the  most  effect  on  the  response  of  the  simulation.  Any  increase  from  unity 
resulted  in  oscillations  in  the  related  control  surfaces.  Since  the  maximum  angles 
that  the  thrust  actuators  can  move  is  20  degrees  and  the  maximum  angle  the  control 
surface  can  move  is  25  degrees,  the  thrust  actuators  vectoring  saturate  before  the 
aerodynamic  surfaces.  Because  of  this,  gains  on  thrust  vectoring  are  lowered  below 
the  level  that  causes  control  surface  oscillations  to  help  prevent  these  saturations. 
Another  indicator  for  the  thrust  gains  is  the  amount  of  Note,  as  the  gain  on  the 
yaw  channel  of  thrust  vectoring  is  reduced  from  unity  the  amount  of  ^  decreases 
until  a  particular  gain  is  reached  and  then  0  increases. 

4-10  Maneuvers 

The  maneuvers  are  produced  by  trial  and  error.  An  input  vector  that  repre¬ 
sents  the  command  for  the  desired  output  is  formulated  using  Matrixx  and  then  the 
simulation  is  executed  using  this  input  vector.  The  resulting  output  is  plotted  versus 
time  and  the  input  is  adjusted  until  the  desired  output  is  produced.  To  represent 
the  specifications  on  the  pitch-up/pitch-down,  four  maneuvers  are  chosen  and  imple¬ 
mented  as  shown  in  Figure  3.6  through  3.8.  To  represent  all  of  the  conditions  m  this 
envelope,  each  pitch  down  rate  must  be  accomplished  after  each  specific  pitch-up 
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maneuver.  This  gives  sixteen  responses.  When  the  nose  of  the  aircraft  is  not  pitched 
up  very  high  the  fast  pitch  down  causes  problems  with  the  system  identification 
which  is  discussed  later,  and  so  these  pitch  down  maneuvers  are  not  implemented. 
Three  roll  maneuvers  are  chosen  to  represent  the  roll  specifications  as  shown  in  Fig¬ 
ure  3.9  through  3.11.  Each  of  these  roll  maneuvers  must  be  accomplished  from  the 
pitch  attitude  developed  from  one  of  the  four  pitch-up  maneuvers  so  there  are  twelve 
maneuvers  associated  with  roll.  All  together  there  are  twenty-two  maneuvers  to  rep¬ 
resent  the  specifications.  At  this  point  the  input/output  vector  time  histories  are 
stored  in  file  to  be  used  in  the  system  identification  process. 


4.11  Initial  System  Identification 

As  stated  before  the  equivalent  plant  is  found  through  the  use  of  maxlike  in 
MatriXx  .  The  first  model  used  to  define  the  equivalent  plant  is  shown  in  Equations 
(4.1)  through  (4.3). 
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Pij  =  +  Pb.,S^  +  Pci,S^  +  Pdi,S  -f  Pe.j 

di  =  P/.S®  -f-  Pg^S^  -f-  Phi^^  +  Pi,S^  +  Pj,S  -f  Pki 


This  model  is  a  matrix  of  transfer  functions  each  having  fourth-order  over  fifth-order, 
numerator  over  denominator,  with  a  common  denominator  for  each  row.  This  model 
proves  to  be  unsatisfactory  because  of  non-minimum-phase  (NMP  )  plants,  which  is 
explained  later  in  this  chapter.  From  earlier  testing  the  initial  parameter  vector  is 
chosen  to  be  a  vector  of  zeros.  Results  from  this  setup  are  not  good,  the  equivalent 
plant  does  not  reproduce  the  original  desired  output  with  the  corresponding  input. 


4.12  Time  Segments 

After  trying  different  initial  parameters,  with  no  success,  a  smaller  time  seg¬ 
ment  is  used  to  see  if  smaller  individual  features  can  be  modeled  and  the  resulting 
equivalent  plant  used  as  the  initial  parameter  vector  for  a  larger  segment  of  time. 
The  results  are  mixed,  the  smaller  segment  of  time  produces  an  acceptable  model, 
one  that  follows  the  major  features  of  the  response,  but  when  this  plant  is  used  as 
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an  initial  vector  the  equivalent  plant  is  not  acceptable.  At  this  point  it  is  decided 
that  since  there  are  intervals  of  time  where  the  dynamics  of  the  maneuvers  are  close 
to  steady  state  which  results  in  very  small  derivatives  in  these  intervals,  the  ma¬ 
neuvers  can  be  broken  into  time  segments  and  used  in  the  nonlinear  QFT  process 
as  individual  plants.  This  time  division  is  made  at  the  three  second  point  in  all 
maneuvers.  The  first  three  seconds  of  all  of  the  maneuvers  is  a  pitch  up.  For  the 
nonlinear  QFT  process  the  four  different  pitch-up  maneuvers  result  in  four  plants. 
The  ten  different  pitch-down  maneuvers  are  separate  plants  as  are  the  twelve  roll 
maneuvers.  Subtracting  the  six  plants  not  modeled  due  to  the  pitch-down  problem 
gives  a  total  of  twenty  equivalent  plants  to  model. 

4.13  Modeling  Difficulties 

The  model  chosen  to  represent  the  nonlinear  aircraft  equations  is  a  matrix  of 
transfer  functions  that  are  each  fourth-order  over  fifth-order  as  shown  in  Equations 
(4.1)  and  (4.2).  Modeling  the  remaining  pitch-up/pitch-down  maneuvers  is  accom¬ 
plished  with  no  difficulties.  The  roll  maneuver  presents  some  difficulty,  probably 
because  of  the  coupled  nature  of  the  maneuver.  To  try  to  get  a  good  transfer  func¬ 
tion  fit,  different  initial  parameters  are  chosen.  To  choose  the  different  parameters, 
knowledge  of  the  physics  of  the  aircraft  and  an  insight  in  the  operation  of  transfer 
functions  are  used  to  separate  the  transfer  function  matrix  into  three  sub-matricies. 
The  idea  is  to  find  the  parameters  for  these  sub-models  and  then  combine  the  re¬ 
sults  into  an  initial  parameter  vector  for  the  whole  physical  model.  This  approach 
produces  acceptable  transfer  function  fits. 

4.14  Non-Minimum-Phase  Plants 

After  the  equivalent  plants  have  been  generated,  they  are  examined  by  factoring 
the  polynomials  into  roots.  At  this  point  the  equivalent  plants  are  found  to  be  NMP. 
These  zeros  are  in  the  range  of  1  to  10,  this  means  that  the  QFT  design  technique 
cannot  be  used  for  these  transfer  functions.  According  to  Horowitz,  the  only  control 
design  method  that  has  been  shown  capable  of  a  design  with  this  relative  magnitude 
of  right-half-plane  (RHP  )  zeros  is  the  singular-g  method  [9].  The  singular-g  method 
was  used  by  Walke  in  his  thesis  and  essentially  uses  matrix  manipulation,  a  singular 
compensation  matrix,  and  trial  and  error  to  design  a  compensator.  To  try  to  prevent 
NMP  zeros,  the  initial  parameters  are  manipulated  to  try  to  keep  all  of  the  zeros 
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negative,  this  does  not  work.  After  much  work  in  manipulating  the  parameters,  a 
new  model  is  developed  that  constrains  the  equivalent  plants  to  be  MP.  This  model, 
shown  in  Equations  (4.4)  and  (4.5)  ,  uses  MatriXx  functions  such  as  the  absolute 
value,  abs,  to  insure  no  RHP  zeros  are  allowed  in  the  equivalent  plants. 
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4-15  Non-Minimum-Phase  Effects 

There  are  parts  of  the  maneuvers  that  can  not  be  modeled  under  the  constraint 
that  a  MP  model  be  achieved.  By  removing  the  /3  inner  loop  feedback  it  is  possible 
to  model  the  roll  maneuvers.  It  appears  that  the  ^  feedback  causes  the  plants  to 
be  strongly  NMP.  The  /3  feedback  is  removed,  as  shown  in  Figure  4.1,  and  all  of 
the  maneuvers  are  simulated  to  collect  the  input/output  tin.e  histories.  Six  of  the 
maneuvers  can  not  be  modeled  by  this  method.  These  maneuvers  represent  the  fast 
pitch-down  after  a  relatively  slow  pitch-up.  It  appears  that  these  maneuvers  can  not 
be  modeled  as  MP  and  the  decision  is  made  to  exclude  them  from  the  flight  envelope. 
This  decision  is  justifiable  because  these  maneuvers  have  the  aircraft  pitching  down 
well  past  the  attitude  intended  by  the  concept  of  the  maneuvers. 


4-16  Non-Minimum-Phase  q’s 

After  finding  the  equivalent  MP  plants  for  the  maneuvers,  it  is  necessary  to 
invert  the  P  matrices  to  obtain  the  P*  matrices.  Upon  completion  of  this  inversion, 
the  individual  qij(o)  transfer  functions  are  found  to  be  NMP.  The  solution  to  this 
problem  is  to  bypass  the  direct  development  of  the  P  matrices  and  use  the  system 
identification  routines  to  directly  identify  the  q,jS  under  the  constraint  of  achieving 
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MP  qijS,  see  Equations  (4.6)  and  (4.7). 
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where, 


,  _  J_  _  Pg.j-S  +  PbijS‘^  +  PcjjS  +  Pd,jS  +  fe.j -S  + 

^g.j('5+  Ph,j  )(5+  Pi^j  )(s+  Pj^j  )(s+  Pk,j  )(s  +  100) 


There  are  benefits  derived  from  this  solution.  One  is  that  no  inversion  of  the  three  by 
three  matrix  of  transfer  functions  is  needed,  easing  the  numerical  burden  on  MatriXx  . 
Another  is  lifting  the  constraint  that  the  rows  in  the  transfer  function  matrices  had  to 
have  a  common  denominator.  Modeling  the  equivalent  plant  as  a  matrix  of  transfer 
functions  with  individual  denominators  produces  better  fits  and  the  time  required 
for  system  identification  is  reduced  dramatically.  Plots  of  the  input,  equivalent  plant 
input,  and  outputs  are  shown  in  Appendix  2.  All  of  the  plants  are  modeled  and  the 
transfer  functions  are  shown  in  Appendix  3. 


4-17  Conclusion 

This  chapter  covers  the  development  of  equivalent  plants  that  are  used  in  the 
nonlinear  QFT  method  for  FCS  design.  First  system  identification  is  discussed 
with  a  concentration  on  batch  identification  techniques.  This  is  followed  by  by  a 
discussion  of  the  methods  used  to  obtain  the  input/output  time  histories  needed  for 
system  identification.  Problems  with  identifying  the  plants  is  discussed  including 
NMP  plants,  long  time  segments,  and  NMP  qiS.  In  the  next  chapter,  the  method  of 
finding  the  effective  plants  and  Matrixx  generated  boundaries  are  discussed. 


4-9 


V.  QFT  FCS  Design 


5. 1  Introduction 

This  chapter  details  the  development  of  Qii{k)  and  the  design  of  boundaries  that 
dictate  the  shape  of  The  first  section  of  the  chapter  explains  an  experiment 

showing  that  Matrixx  does  not  have  a  robust  method  of  accommodating  high  order 
rational  transfer  functions,  which  has  an  observable  impact  on  the  development  of 
Qij{k)-  Next  is  a  discussion  of  Qij{k)  and  how,  in  some  cases,  boundaries  can  be 
developed  to  reduce  the  number  of  its  large  RHP  poles.  The  final  topic  of  this 
chapter  is  the  development  of  the  stability  boundaries  and  the  boundaries  used  to 
reduce  large  RHP  poles  which,  in  this  thesis,  are  called  P-boundaries  . 

5.2  Matrixx  Algorithm  Robustness 

During  the  development  of  Qij(k)i  plants  are  developed  that  are  as  high  as  54*^ 
order.  These  high  order  transfer  functions  are  probably  due  to  inaccuracies  which  are 
a  result  of  the  MatriXx  algorithms  that  are  not  robust  enough  to  accommodate  high 
order  plants.  A  simple  "'xperiment  is  performed  to  show  the  limits  of  the  MatriXx 
algorithm  robustness  [20].  Ten  eigenvalues  are  taken  whose  values  range  from  -1  to 
-10.  These  eigenvalues  are  used  to  form  the  polynomial  B,  as  shown  in  Equations 
(5.1)  and  (5.1). 

A  =  (A  -1-  1)(A  +  2)(A  -h  3)(A  -h  4)(A  -h  5)(A  -f  6)(A  +  7)(A  -h  8)(A  +  9)(A  -b  10)  (5.1) 

B  =  polyi^A)  =  A  ^  "b  flgA^  -b  (igX^  +  djXJ  -b  o,q)^  -b  OsA^  -b  -b  cl^X?  -b  (I2X?  -b  OjA^  -b  clq 

(5.2) 

This  is  accomplished  in  Matrixx  by  entering  the  vector,  A‘=[-l  -2  -3  -4  -5  -6  -7  -8 
-9  -10],  whose  elements  reresent  the  eigenvalues.  The  polynomial  is  then  formed  by 
using  the  Matrixx  command,  B=poly(a).  The  location  of  the  roots  of  the  resulting 
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polynomial  are  plotted  in  Figure  5.1.  As  seen  in  the  figure,  some  warping  has  taken 
place,  that  is,  the  given  real  eigenvalues  have  migrated  into  the  imaginary  region.  By 
changing  the  numbers  of  roots  it  is  found  that  plants  with  8*^  order  polynomials  do 
not  exhibit  the  warping  associated  with  manipulations  of  higher  order  plants.  As  the 
order  of  plants  increase  the  magnitude  of  the  warping  seems  to  increase,  see  Figure 
5.2  for  a  54*^  order  example.  As  in  the  first  example,  the  second  example  starts  with 
a  54*^  order  vector  of  all  real  eigenvalues,  this  time  the  A^’s  are  arranged  from  -1  to 
-54  at  increments  of  -1  along  the  real  axis.  When  the  eigenvalues  are  transformed 
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Figure  5.1.  Matrixx  Robustness  10*^  Order  Real  Roots 

into  polynomial  form  as  previously  discussed,  the  roots  of  B,  a.s  shown  in  Figure  5.2, 
result.  The  implication  of  this  experiment  is:  manipulations,  such  as  convolution  of 
transfer  functions  result  in  transfer  functions  where  t’^ere  are  not  many  pole/zero 
cancellations  due  to  warping,  whereas  in  reality  there  can  be  common  poles  and 
zeros  in  these  plants  that  should  cancel.  This  compounds  the  problem  because  plant 
manipulations  result  in  abnormally  high  order  plants  and  any  subsequent  manipula¬ 
tions  of  these  plants  result  in  even  more  severe  warping  of  poles  and  zeros.  Matrix* 
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Figure  5.2.  Matrixx  Robustness  54‘^  Order  Real  Roots 

is  used  for  all  the  manipulations  in  this  thesis  because  this  problem  did  not  show  up 
until  late  in  the  design  process.  Subsequent  theses  will  need  to  consider  this  problem 
early  in  the  process. 


5-3  Qij{k)  Development 


According  to  Horowitz,  to  ensure  stability  for  the  MIMO  problem,  Qij(k)  must 
be  developed  [9]. Since  this  thesis  only  uses  stability  bounds,  cross-coupling  in  the 
plant  matrix  is  not  accounted  for  in  disturbance  boundaries,  therefore  the  only  factor 
that  accounts  for  these  cross-coupling  terms  is  the  use  of  Qij{k)  in  the  design  process. 
The  equations  used  to  develop  repeated  from  Chapter  2,  are  Equations  (5.3) 

and  (5.4). 


_  go(fc-i)(l  +  h(k-i)) 
_  <Iik{k-i)qk3{k-i) 


(5.3) 

(5.4) 
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These  equations  are  programmed  in  Matrix*  to  develop  Qij(k)-  As  explained  in 
Chapter  6,  the  overall  FCS  is  divided  into  FCSs  that  are  scheduled  in  the  nonlinear 
simulation.  Each  FCS  is  designed  using  the  same  method  as  the  others,  therefore 
this  description  concentrates  on  the  design  of  the  first  loop  and  only  give  details  of 
the  design  of  the  other  loops  when  important.  After  /o,  has  been  designed,  Equations 
(5.3)  and  (5.4)  are  used  to  find  922(1)  and  933(1)  for  each  plant  to  be  controlled  by 
this  FCS,  these  are  the  sets  ^22(1)  and  Q33(i).  To  find  033(2),  the  sets  of  plants: 
032(1))  023(1))  and  022(1)  also  need  to  be  calculated.  With  these  plants  calculated 
the  set  033(2)  is  calculated.  When  a  comparison  is  made  of  the  Bode  plots  of  022(i) 
and  033(2)  against  022(o)  and  033(o))  it  is  apparent  that  there  are  large  RHP  poles 
in  033(2)  that  do  not  appear  in  033(o).  This  is  a  problem  because  these  large  RHP 
poles  cause  uncertainty,  in  phase  angle  of  the  frequency  response  of  033(2))  to  be 
greater  than  300°  until  almost  10,000  ^  .  This  causes  the  phase  margin  frequency 
to  be  at  least  10,000  ^  .  This  situation  is  also  found  in  the  second  and  third  loops 
of  the  second  FCS.  If  certain  conditions  are  met,  these  large  RHP  poles,  induced  by 
the  effective  plant  process,  can  be  prevented.  Since  the  third  loop  is  modified  by  the 
two  previous  loops,  if  these  loops  can  be  shaped  in  a  manner  that  make  Equations 
(5.5)  and  (5.6)  true,  or  if  Equation  (5.5)  is  true  and  Equation  (5.7)  is  true  for  all  of 
the  plants,  then  the  large  RHP  poles  can  be  prevented  [9].  In  these  equations  is 
some  positive  number. 


or 


1  -  7.i(fc+i)  +  ^(fe) 


(jis  ^  <  0 


(5.5) 

(6.6) 


(5.7) 


Equation  (5.5)  is  used  to  prevent  excursions  across  the  ju;  axis,  which  result  in  large 
RHP  poles.  Equation  (5.6)  is  a  condition  that  Equation  (5.5)  must  satisfy  unless  all 
of  the  plants  satisfy  Equation  (5.7).  These  are  necessary  conditions  for  even  moderate 
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feedback  benefit.  These  equations  axe  implemented  through  the  use  of  T-boundaries, 
which  is  discussed  later.  When  the  sets  ^33(1)  and  Q33(2)  are  examined  it  is  found 
that  large  RHP  poles  are  induced  by  the  set  of  /i(o)S  (-Ci(o))  in  the  denominator 
of  Equation  (5.3)  and  no  large  RHP  are  induced  by  jC2(i).  For  the  second  PCS, 
RHP  poles  are  induced  by  both  £1(0)  and  £2(1)-  When  Equation  (5.6)  is  evaluated 
based  on  Q33(i)  and  w=10®  ^  there  are  two  ^33(1)  that  fail  the  condition.  The 
plants  corresponding  to  these  two  533(1)  can  either  be  eliminated  or  a  separate  ECS 
design  can  be  made  for  them.  The  decision  is  made  to  eliminate  these  two  plants  and 
continue  with  the  design  process.  This  decision  is  made  primarily  based  on  the  added 
time  to  design  another  ECS  and  that  these  plants  are  not  plants  on  the  boundaries 
of  the  maneuver  envelopes.  For  the  second  PCS,  it  is  found  that  6  533(1)  out  of  12 
violate  Equation  (5.6).  Based  on  the  Bode  plots  of  these  plants,  it  is  decided  to 
split  the  third  channel  of  the  second  ECS  into  three  smaller  FCSs,  this  is  discussed 
further  in  Chapter  6.  By  reshaping  the  /oi(o)’s,  based  on  boundaries  developed  from 
Equation  (5.5)  and  trial  and  error,  all  of  the  large  RHP  poles  induced  in  ^33(1)  are 
prevented.  After  shaping  the  necessary  finding  Qij(k)  is  a  matter  of  running  the 
MatriXx  programs  to  do  the  calculations.  Because  of  the  problem  of  robustness  of  the 
Matrix*  algorithms,  two  different  methods  of  transfer  function  reduction  are  used, 
with  very  limited  results.  One  of  the  methods  is  a  pole-zero  cancellations  program 
written  by  Miller,  [18],  that  uses  Matrix*  functions  to  separate  the  transfer  function 
polynomials  into  roots  and  then  compares  the  roots  for  cancellation.  The  second 
method  is  a  Matrix*  function  called  ’minimal’  [16].  Neither  method  overcomes  the 
warping  problem  and  this  results  in  very  high  order  plants  which  are  suspect. 

5.4  lo  Boundaries 

Because  of  the  specification  that  the  phase  margin  frequency  be  no  more  than 
30  ^  ,  Horowitz  points  out  that  the  stability  boundaries  will  tend  to  be  the  dominate 
boundaries  and  it  is  decided  to  base  the  design  on  stability  boundaries  alone.  To 


5-5 


develop  the  stability  boundaries,  templates  representing  the  magnitude  and  phase 
uncertainty  of  the  plants  are  developed  using  Matrixx  and  plotted  on  paper.  The 
idea  is  to  move  these  templates  around  the  3  dB  M-contour  on  a  Nichols  chart, 
while  tracking  a  nominal  plant  {qoij(k))‘  The  path  that  qoij(k)  traces  on  the  chart  is 
the  stability  boundary.  The  3dB  M-contour  is  chosen  because  it  represents  a  phase 
margin  of  35  degrees.  Instead  of  moving  the  template  on  the  Nichols  chart,  Lacey 
developed  a  method  in  which  the  Nichols  chart  with  the  3dB  M-contour  is  plotted  on 
clear  acetate  and  the  template  is  held  stationary  while  the  Nichols  chart  is  moved  fl3]. 
After  developing  stability  boundaries  for  the  diagonal  plants,  the  problem  with  large 
RHP  poles  is  encountered.  These  RHP  poles  make  it  necessary  to  develop  a  MatriXx 
program  that  plots  the  P-boundaries  and  redetermines  the  stability  boundaries. 
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Figure  5.3.  Example  Template,  u;=100  — 


5.5  MatriXx  Generated  Stability  Boundaries 

Horowitz  devised  a  method  in  which  Equation  (5.8)  can  be  used  in  a  computer 
program  to  develop  stability  boundaries,  where  Mm  is  the  inverse  log-magnitude  of 
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the  desired  M-contour,  in  this  case  M^.  =  1.4125,  and  £,(*;) 
shaped. 


1  +  A(fc) 


<  Mm 


is  the  set  of  loops  to  be 


(5.8) 


This  method  allows  boundaries  on  the  magnitude  of  the  nominal  loop  {lonk)) 
found  as  a  function  of  its  phase.  To  develop  this  program,  the  following  equation  is 
formed  that  relates  £,  Qij(k)i  and  the  nominal  plant  {qoij{k))' 


Qim 


Equation  (5.8)  is  rearranged  resulting  in  Equation  (5.10). 


(5.9) 


Mm 


1  +  ^t(A:) 


> 


(5.10) 


Both  sides  of  Equation  (5.10)  are  squared,  and  Equations  (5.11)  and  (5.12)  are 
substituted  in,  where  m  is  the  magnitude  of  hnk-t)  and  6  is  the  phase.  This  results 
in  Equation  (5.13),  where  A,  B,  and  C  are  defined  by  Equations  (5.14)  through 
(5.16). 


=  m(cos(6>)  +  jsin(0)) 

(5.11) 

Qij{k)  , 

=  =  f^2  +  J  ^2 

(5.12) 

Am^  +  Bm  +  C  >  0 

(5.13) 

Ml  +  fil) 

(5.14) 

B  =  [2M^  ifii  cos{e)  -  H2  sin(0))  -  +  fi])] 

(5.15) 

C=Mm 

(5.16) 

By  solving  this  quadratic  equation  in  m  for  a  specific  6,  w,  %(fc),  and  Qc 

,  two  roots 

are  found  that  correspond  to  the  magnitude  of  the  boundaries  on  . 

This  equation 
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must  be  solved  for  each  qij(k)  £  Qij{k)i  for  each  0,  for  the  boundaries  at  one  value  of 
uj.  The  upper  stability  boundary  magnitude  for  that  0  is  the  maximum  value  from 
the  set  of  the  maximum  values  of  the  two  roots  for  all  Qij(k)-  The  lower  boundary  is 
the  minimum  value  from  the  set  of  minimum  values  for  all  Qij(k)-  Example  plots  of 
boundaries  are  shown  in  Figures  5.4  through  5.4.  Comparing  the  Matrix*  generated 
boundaries  and  the  hand  drawn  boundaries  shows  good  correlation.  A  good  test  of 
these  boundaries  is  the  plot  of  the  set,  on  a  Nichols  chart  as  in  the  example  in 
Figure  5.8. 


Figure  5.4.  Example  Stability  Bound,  for  at  a;=l  through  10  ^ 
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Figure  5.7.  Example  Stability  Bound  £©2(1)  at  u;=10  ^ 
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Figure  5.8.  Example  Plot  of  £1(0)  on  Nichols  Chart 
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5.6  Matrixx  Generated  F-Boundaries 


Because  of  the  necessity  to  prevent  the  development  of  large  RHP  poles  in 
the  effective  plant  process,  Horowitz  developed  the  following  method  to  mechanize 
F-boundaries.  The  development  is  directly  analogous  to  the  development  of  the 
stability  boundaries  and  is  based  on  satisfying  Equation  (5.5),  which  is  repeated  in 
Equation  (5.17),  in  this  case,  the  choice  of  —  0.5  is  made. 


1  -  7ij(fc+i)  +  h(k) 


>^P 


(5.17) 


As  in  the  case  of  the  stability  bounds,  the  object  is  to  find  boundaries  on  the  mag¬ 
nitude  of  35  a  function  of  its  phase.  Both  sides  of  Equation  (5.17)  are  squared 
and  Equations  (5.11),  (5.12),  and  (5.18)  are  substituted  in  the  resulting  Equation 
(5.19)  where  A,  B,  and  C  are  defined  by  Equations  (5.20)  through  (5.22). 


7.i(fc+i)  =  7i  +  i72  (5.18) 

-I- Hm  +  C  >  0  (5.19) 

^=(7?  +  7|)  (5.20) 

B  =  2[(1  -  7i)  (^1  cos{$)  -  H2sm{0))  -  72  {ni  sin(0)  +  /X2cos(0))]  (5.21) 

C  =  (1-7i)'  +  72-^'  (5.22) 


The  F-boundaries  are  then  generated  in  the  same  manner  as  the  stability  boundaries. 
Example  plots  of  F-boundaries  are  shown  in  Figures  5.9  through  5.12.  A  test  of  these 
boundaries  is  if  large  RHP  poles  are  generated  after  the  use  of  F-boundaries.  In  most 
cases  no  large  RHP  poles  are  generated,  but  there  are  two  cases  where  large  RHP 
poles  are  generated.  In  these  two  cases  the  are  manipulated  by  trial  and  error 
until  no  large  RHP  poles  are  generated.  It  is  unclear  in  both  of  these  cases  where  the 
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Figure  5.9.  Example  F-Boundaries  for  .)  at  uj=1  ^  ,  Based  on  033(0) 

violations  of  the  F-boundaries  occurs,  but  in  all  the  others  cases  where  large  RHP 
poles  are  gei'.erated,  clear  violations  of  the  F-boundaries  are  present. 

5.7  Matrixx  Generateri  Boundaries:  General  Considerations 

When  using  these  algorithms  to  generate  stability  boundaries  and  F-boundaries, 
there  are  many  factors  to  consider.  In  this  Matrix*  program  the  boundaries  are  con¬ 
strained  to  start  and  stop  at  the  0  dB  axis.  This  causes  the  boundaries  to  be  exag¬ 
gerated  at  the  ends.  When  there  are  fast  changes  in  the  boundaries,  the  resolution 
along  the  phase  axis  is  a  major  factor  in  how  exaggerated  the  boundaries  become. 
When  there  is  large  uncertainty  in  the  phase  of  the  %(*),  see  Figure  5.13,  separate 
regions  of  boundaries  are  formed.  In  reality  these  regions  should  be  connected  with 
a  continuous  line,  and  care  must  be  taken  to  identify  the  regions  and  shape 
completely  around  the  right  side  of  the  boundaries.  Boundaries  of  plants  that  have 
large  phase  uncertainty  can  cause  the  width  of  the  boundary  to  be  masked  by  the 
repetitive  nature  of  the  Nichols  chart.  That  is,  for  phase  uncertainty  more  than  360° 
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Figure  5.10.  Example  F-Boundariesfor  at  u;=l  ^  ,  Based  on  622(1) 

the  boundaries  on  the  far  right  edge  of  the  boundary  overlap  onto  the  region  be¬ 
tween  -180®  and  -0°.  These  overlapped  boundaries  give  a  false  indication  of  the  right 
edge  of  the  boundary.  This  matter  is  complicated  when  there  are  large  separations 
between  the  phase  of  the  plants.  The  area  between  the  separations  do  not  result  in 
boundaries.  This  cause  boundaries  with  numerous  non-contiguous  regions  as  shown 
in  Figures  5.14  through,  5.21.  The  min’mum  frequency  at  which  a  loop  can  cross 
the  0  dB  axis  in  these  figures  is  500  ^  .  As  shown,  the  complete  right  edge  of  the 
boundary  emerges  from  the  boundaries  around  the  prohibited  region  that  is  360°  to 
the  right  of  the  original.  This  perception  problem  can  be  prevented  by  algorithmi- 
caiy  drawing  a  closed  contour  around  the  mathematical  template,  as  suggested  by 
Sating  [20]. 
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ire  5.11.  Example  F-Boundaries  for  £02(1)  a;=100  ^  ,  Based  on  Q22(i) 


Figure  5.12.  Example  F-Boundaries  for  £02(1)  at  u;=50  ^  ,  Based  on  ^22(1) 
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Figure  5.15.  Large  Q  Uncertainty  Example  cont’d,  Stability  Boundaries  for  £, 
at  u;=50  ^ 


Figure  5.16.  Large  Q  Uncertainty  Example  cont’d,  Stability  Boundaries  for 
at  cj=70  — 
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Figure  5.17.  Large  Q  Uncertainty  Example  cont’d,  Stability  Boundaries  for  >^03(2) 
at  a;=90  ^ 
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Figure  5.18.  Large  6  Uncertainty  Example  cont’d,  Stability  Boundaries  for  £03(2) 
at  a;=200  ^ 
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Figure  5.19.  Large  0  Uncertainty  Example  cont’d,  Stability  Boundaries  for  C 
at  a;=500  ^ 
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Figure  5.20.  Large  Q  Uncertainty  Example  cont’d,  Stability  Boundaries  for  £ 
at  u;=700  ^ 
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Figure  5.21.  Large  Q  Uncertainty  Example  cont’d,  Stability  Boundaries  for  C 
at  u;=1000  ^ 


5.8  Conclusion 


This  chapter  covers  development  of  the  effective  plants  and  generation  of 
boundaries  on  First,  limitations  of  the  Matrixx  algorithms  for  manipulating 

plants  are  discussed  and  a  basic  experiment  to  show  this  is  detailed.  The  develop¬ 
ment  of  Qij(k)  and  some  of  the  problem  encountered  is  discussed.  These  problems 
include  the  development  of  large  RHP  poles  by  the  effective  plant  process,  which  is 
solved,  in  most  cases,  by  P-boundaries.  Finally,  a  discussion  of  MatriXx  generated 
boundaries  is  presented.  Horowitz’s  method  for  generating  stability  bounds  is  dis¬ 
cussed  first  followed  by  his  method  for  generating  F-boundaries.  The  last  sections 
outline  some  of  the  cautions  in  using  the  Matrixx  generated  boundaries.  In  Chapter 
6  loop  shaping  and  filter  design  are  discussed. 
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VI.  Loop  Shaping  and  Filter  Design 


6.1  Introduction 

This  chapter  details  the  design  of  the  loop  transmission  and  the  prefilter.  First, 
general  considerations  and  the  reason  for  controller  scheduling  are  presented.  Next, 
design  of  the  loop  transmissions  for  all  FCSs  is  described,  including  many  problems 
with  RHP  poles.  Finally  there  is  a  description  of  the  prefilters  designed  to  cause  the 
aircraft  to  respond  in  a  specified  manner  to  specified  control  inputs. 

6.2  General  Considerations 

After  developing  the  Qij(o),  from  the  system  identification  algorithms,  a  Bode 
plot  is  obtained  for  all  ^,,(0)  €  Qu(o)  i=l,2,3.  Figures  6.1  through  6.3.  On  inspection 


Omega(rad\aec) 


Figure  6.1.  Bode  Plots  for  All  qu[o) 


the  Bode  plots  reveal  that  the  magnitudes  have  an  extreme  amount  of  uncertainty, 
but  the  qa^o)  are  grouped  in  two  groups.  This  indicates  that  scheduling  of  gain  is 
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Figure  6.2.  Bode  Plots  for  All  922(0) 

needed  for  the  FCS  design.  When  the  two  groups  are  separated  and  plotted  on 
Bode  plots,  Figures  6.4  through  6.9,  it  is  evident  that  the  uncertainty  in  phase  is 
too  large  for  just  gain  scheduling  and  separate  controllers  must  be  designed  for  the 
two  groups.  The  decision  is  made  to  develop  two  FCS’s  and  schedule  the  FCS  based 
on  the  maneuver.  Upon  inspection  it  is  found  that  the  maneuvers  associated  with 
pitch-  up/pitch-down  can  be  grouped  together  and  the  maneuvers  representing  the 
velocity  vector  roll  can  be  grouped  together. 

6.S  Loop  Shaping 

The  first  forward  loop  is  shaped  using  hand  generated  stability  boundaries.  A 
Matrix*  program  is  used  to  ease  the  task  of  loop  shaping.  Matrix*  has  an  interactive 
Nichols  chart,  which  means  an  w  of  interest  on  can  be  identified  by  a  cursor 
that  moves  on  Because  of  this,  one  effective  method  of  loop  shaping  is  to  use 

a  Matrix*  program  to  build  a  based  on  poles  and  zeros  and  then  plot  it  on  the 
Nichols  chart.  In  this  manner  can  be  formed  by  adding,  deleting,  or  moving 
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Figure  6.3.  Bode  Plots  for  All  933(0) 


poles  and  zeros  and  then  checking  the  results  on  the  Nichols  chart.  To  do  this  loop 
shaping,  using  the  Nichols  chart,  basic  pole/zero  relationships  must  be  used.  Some 
of  the  more  important  relationships  are: 


•  adding  a  pole  will  pull  to  the  left  and  causes  the  magnitude  to  decrease. 

•  adding  a  zero  will  pull  to  the  right  and  causes  the  magnitude  to  increase. 

•  by  adding  a  pole/zero  pair  can  be  moved  to  the  right  or  to  the  left  without 
changing  the  magnitude  at  the  a;  of  interest.  This  can  be  accomplished  by 
picking  the  position  of  the  pole/zero  and  then  placing  the  zero/pole  an  equal 
distance,  on  a  log  scale,  from  the  of  interest.  The  direction  that  the  loop  is 
moved,  is  left  or  right  depending  on  whether  the  pole  or  the  zero  is  of  lesser 
magnitude 

•  A  complex  pair  of  zeros  or  poles  can  be  used  where  a  large  change  in  phase 
is  needed  quickly  without  a  corresponding  quick  change  in  magnitude.  The 
change  in  phase  and  magnitude  both  depend  on  the  damping  factor  thus  the 
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Figure  6.4.  Bode  Plots  for  FCS#1  9n(o) 


use  of  a  chart,  where  complex  roots  are  plotted  with  changes  in  damping  factor, 
is  helpful. 

•  When  using  a  complex  pair  of  poles  to  swing  lo,^,  iround  the  bottom  of  the 
stability  boundaries,  the  damping  ratio  has  a  m.  ''ffect  on  the  shape  of  the 
corner.  A  large  damping  ratio  causes  the  con  be  smooth  and  displaced 
to  the  left  when  compared  to  the  corner  produced  by  a  low  damping  factor.  If 
the  corner  needs  to  be  moved  right  a  very  low  damping  factor  can  be  used.  A 
small  damping  factor  (.4-. 6)  is  useful  in  turning  the  loop  below  the  stability 
bound. 


These  are  only  the  more  important  relationships  that  are  helpful  in  loop  shaping 
using  the  Nichols  chart.  In  addition  to  these  relationships  one  must  know  when  these 
effects  should  take  effect.  These  relationships  can  be  derived  from  Bode  plot  rules. 
[5]  There  is  a  problem  when  using  the  frequency  response  algorithms  of  Matrix*  , 
and  other  CAD  packages,  for  functions  that  contain  RHP  poles  and/or  zeros,  '’^hen 
RHP  factors  are  present  Matrix*  shifts  the  phase  plots  by  some  value  of  phase.  The 
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Figure  6.5.  Bode  Plots  for  FCS^l  922(0) 


phase  can  be  corrected  by  shifting  the  frequency  data,  in  phase,  by  a  factor  obtained 
by  computing  the  difference  between  the  number  of  poles  and  zeros,  multiplying 
this  number  by  90°,  and  subtracting  the  final  phase  of  the  frequency  data  from  this 
number.  The  result  is  a  phase  correction  factor  that  must  be  added  to  the  phase 
plot. 

6.4  FCS:^1  Loop  Design 

For  of  the  first  FCS,  a  9o,,(o)  is  chosen  to  be  the  lowest  magnitude  for  most 
of  the  frequencies  of  interest.  This  choice  is  made  to  cause  the  boundaries  to  be  below 
the  3  dB  M-contour.  The  drawback  of  this  choice  is  it  makes  it  necessary  that 
be  well  below  the  30  ^  specification  for  phase  margin  frequency(u;^  )  so  that  all  the 
other  /i(o)  G  C\  satisfy  the  <  30^  specification.  With  this  choice,  the  stability 
boundaries  are  formed  by  the  method  of  templates.  Since  9o,i(o)  contains  two  RHP 
poles,  they  are  included  in  lo^o)■  From  this  point,  poles  and  zeros  are  added  to  shape 
the  loop.  After  has  been  designed,  the  2^(1)  needed  for  succeeding  loops  are 
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Figure  6.6.  Bode  Plots  for  FCS#1  933(0) 


calculated.  After  shaping  /o^^o)  and  calculating  ^33(2)  it  is  apparent  that  large  RHP 
poles  have  been  added  to  ^33(2)  and  it  is  decided  to  use  the  method  described  in 
Chapter  6  to  prevent  the  formation  of  large  RHP  poles.  By  examining  Q33(i)  and 
Q33(2)?  its  is  apparent  that  the  large  poles  are  a  result  of  /i(o).  Because  two  933(1) 
violate  the  conditions  on  the  F-boundaries,  they  are  removed  from  the  design.  Since 
two  933(1)  represent  pitch-down  manevers  from  a  low  magnitude  pitch-up  maneuver, 
their  removal  from  the  design  changes  the  specified  boundaries.  A  better  method  in 
handeling  this  situation  is  to  design  a  separate  FCS  for  these  two  933(1)  as  is  done 
for  the  third  channel  of  FCS#2.  The  first  loop  is  reshaped  using  Matrix*  generated 
stability  boundaries  and  F-boundaries.  The  results  are  shown  in  Figures  6.10  and 
6.11.  For  this  loop,  the  F-boundaries  did  not  prevent  the  existence  of  all  of  the 
large  RHP  poles  in  ^33(2)  and  it  is  necessary  to  shift  outward  away  from  the 
boundaries  until  the  large  RHP  poles  do  not  appear.  In  this  case  the  boundaries  do 
not  appear  to  be  correct  since  there  are  phases  where  there  is  an  upper  boundary 
below  zero  and  no  apparent  lower  boundary.  As  explained  later,  the  boundaries  are 
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NomintkJ  PJant  Omcget(rad\scc} 

Figure  C.7.  Bode  Plots  for  FCS^2  911(0) 

very  effective  in  reducing  the  RHP  poles  in  FCS  #2.  Once  is  shaped,  effective 
plants  are  calculated,  and  are  shaped.  It  is  not  necessary  to  use  the  F- 
boundaries  when  shaping  ^22(1)  because,  in  this  case  no  large  RHP  poles  are  induced 
by  £2  on  £3.  Figures  6.12  through  6.15  show  and  Appendix  4  contains 

all  the  Gi,  G2,  and  G3  designed  in  this  thesis.  To  check  the  validity  of  the  the 
sets,  £,(fc)  are  plotted  dluiig  with  the  3  dB  M-contor  in  Figures  6.16  through  6.18. 
Note:  in  Figure  6.18,  5  of  the  12  13(2)  can  not  be  plotted  due  to  the  high  order  of 
their  corresponding  933(2)8. 

6.5  FCS#2  Loop  Desiy:n 

The  second  set  of  £,(*:) ’s  is  more  difficult  to  design  because  of  the  large  number 
of  RHP  poles  generated  by  the  effective  plant  process.  After  shaping  it  is  found 
that  large  RHP  poles  are  generated  in  both  Q22(i)  and  633(1).  .\  check  ’.s  made  on 
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Figure  6.8.  Bode  Plots  for  FCS#2  922(0) 
both  for  the  condition  of  Equation  (5.6),  repeated  here  in  Equation  (6.1). 

Jil^  ^  “  7ij(it+t)(iw)  >  0  (6.1) 

The  entire  set  Q22(i)  meets  the  condition,  but  only  6  of  12  of  the  set  Q33(i)  meet 
tfiis  condition.  It  is  decided  to  design  3  separate  FCS  for  the  third  channel, 
is  shaped  in  the  same  manner  as  the  loops  in  FCS#1.  Because  the  F-boundaries 
for  based  on  ^22(1)  m'e  complicated,  is  re-shaped  many  times  before  it  is 
discovered  that  the  r  boundaries,  that  appear  to  have  breaks,  must  be  completely 
encircled  by  Figures  6.19  and  6.20.  Once  is  shaped  to  encircle  the 

r  boundaries,  2-22(1)  c.ii  ulated  and  shaped,  see  Figures  6.21  and  6.22.  The 
shaping  of  is  accomplished  without  the  use  of  F-boundaries  because  large  RHP 
poles  are  not  generated  in  the  933(2)  that  do  not  violate  the  condition  of  Equation 
(6.1).  .\s  stated  before,  233(2)  is  divided  into  three  subsets.  A  separate  FCS  will  is 
desigiK'd  for  ear  h  subset,  d'he  shaping  of  for  each  of  the  subsets  is  accomplished 
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Figure  6.9.  Bode  Plots  for  FCS#2  933(0) 


in  the  same  manner  as  the  previous  loops,  see  Figures  6.23  through  6.25.  The  only 
problem  with  these  loops  is  that  one  of  them  starts  with  6  RHP  poles  which  causes 
it  to  have  a  of  3000  ~  .  £{(*;)  for  FCS#2  are  plotted  in  Figures  6.26  through 
6.30. 
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Figure  6.12.  Loop  Shape  for  FCS#1  /2(i)  with  Stability  Boundaries 
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Figure  6.13.  Loop  Shape  for  FCS#1  l2{i)  with  F-Boundaries 
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Figure  6.14.  Loop  Shape  for  FCS#1  /3(2)  with  Stability  Boundaries 
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Figure  6.15. 


Loop  Shape  for  FCS#1  /3(2)  with  F- Boundaries 


6-12 


In  [dB]  Gain  fdBj 


iMop^HO  .  FC!t  *3 


Figure  6.19.  Loop  Shape  for  FCS#2  ,(o)  with  StabiliL'  Boundaries 


Laopifl  lO  .  FCS 

Figure  6.20.  Loop  Shape  for  FCS#2  /i(o)  with  F-Boundaries 
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Figure  6.21.  Loop  Shape  for  FCS#2  /2(i)  with  Stability  Boundaries 
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Figure  6.22.  Loop  Shape  for  FCS#2  /2(i)  with  F-Boundaries 
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Figure  6.23.  Loop  Shape  for  FCS#21  1^2)  with  Stability  Boundaries 


Phms^  fdc£j 

Loop£33{S)  .  PCS  £33 

Figure  6.24.  Loop  Shape  for  FCS#22  l^2)  with  Stability 
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Figure  6.25.  Loop  Shape  for  FCS#23  1^2)  with  Stability 
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Figure  6.26.  All  FCS^2  /i(o)S 
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Figure  6.30.  All  FCS#22  /3(2)S 

6.6  uj^  Specifications 

Throughout  the  loop  shaping  process,  it  is  difficult,  and  for  many  of  the  loops, 
impossible,  to  meet  the  30  ^  specification  on  uj^  .  The  problem  lies  in  the  number 
of  RHP  poles  involved  in  the  loops.  When  there  are  RHP  poles  present  in  it  is 
necessary  to  start  shaping  with  these  poles.  A  RHP  real  pole  causes  a  phase  of 
—  180°  at  0  ^  .  This  shifts  to  the  left.  For  example,  with  three  RHP  real  poles 
present,  must  translate  at  least  540°  to  the  right  on  the  NC  before  encircling 
the  boundaries.  This  large  translation  causes  the  large  .  For  this  thesis  only  1 
out  of  8  loop  designs  met  the  specified  .  All  maximum  are  shown  in  Table  6.1. 


6.7  Prefilter  Design 

riie  method  used  to  design  the  prefilters  is  the  same  as  outlined  in  Chapter 
2.  I'hrst,  fie  and  B[^  transfer  functions  are  found  that  delineate  the  time  response 
boundaries.  The  time  response  specifications  are  based  on  M1L-STI)-1797A,  ‘Flying 
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Table  6.1.  Maximum  on  All  FCS 


FCS# 

loop  #  (i) 

1 

1 

40 

1 

2 

50 

1 

3 

320 

2 

1 

60 

2 

2 

100 

21 

3 

270 

22 

3 

30 

23 

3 

3000 

Qualities  of  Piloted  Aircraft’.  [17]  Next,  the  frequency  response  of  C^k)  is  plotted 
on  the  same  plot  as  Bu  and  B^.  Finally,  F,-  is  designed  so  that  C^k)  is  contained 
between  Bu  and  Bu-  At  higher  w’s,  Bu  can  be  shifted  lower  by  adding  high  frequency 
poles,  and  therefore  violations  of  Bu  at  high  frequencies  can  be  ignored.  Because 
of  limitations  in  using  MatriXx  ,  some  of  the  prefilters  have  to  be  approximated. 
Figures  6.31  through  6.38  show  the  Bode  plots  of  the  final  designed  CLTFs’.  Some 
of  the  time  responses  can  not  be  plotted  due  to  the  numerical  limitations  of  MatriXx 
and  the  high  order  of  the  qij{k)s.  The  time  responses  that  can  be  plotted  are  shown 
in  Figures  6.39  through  6.43.  Appendix  5  contains  all  the  Fi,  F2,  and  F3  designed 
in  this  thesis. 
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Figure  6.31.  FCS#1  /i(o)S  Filter  Shaping 
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Figure  6.32.  F('S#1  /2(i)S  F'ilter  Shaping 
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Figure  6.39.  FCS#1  /i(o)S  Time  Response 


Figure  6.40.  FCS^l  /2(i)S  Time  Response 
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Figure  6.41.  FCS#2  /i(o)S  Time  Response 
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Figure  6.42.  FCS#22  /3(2)S  Time  Response 
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Figure  6.43.  FCS^23  /3(2)S  Time  Response 
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6.8  Conclusion 


This  chapter  describes  G  and  F  development.  First  the  reasons  for  scheduling 
are  discussed.  Next  the  shaping  of  individual  loops  is  covered.  During  the  loop 
shaping  process,  it  is  found  that  most  of  the  loops  violate  the  specified  because 
of  constriants  on  the  loops.  Finally,  the  design  of  the  prefilters  is  covered.  When 
the  CLTF’s  are  simulated  the  results  are  generally  good.  Chapter  7  described  the 
modeling  problem  that  become  evident  during  the  nonlinear  simulation,  that  prevent 
a  successful  simulation. 
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VII.  Problems  Encountered  with  Nonlinear  Simulation 


7. 1  Introduction 

In  this  chapter,  problems  encountered  during  the  nonlinear  simulation  are  de¬ 
scribed.  The  nonlinear  simulation  is  based  on  Cord‘s  six  degree  of  freedom  nonlinear 
YF-16  simulation  program.  [4]  First,  the  need  for  a  method  of  controller  order  re¬ 
duction  is  outlined.  Next,  a  problem  in  the  modeling  process  is  described.  Finally, 
the  simulation  of  a  flawed  pitch-up/pitch-down  maneuver  is  described.  Note:  the 
pitch-up/velocity  vector  roll  maneuvers  are  not  simulated  due  to  problems  encoun¬ 
tered  with  the  nonlinear  simulation  program  as  explained  at  the  end  of  the  next 
section. 

7.2  Reduced  Order  G,jS 

After  programming  the  designed  FCS  into  the  YF-16  simulation  program,  it 
is  found,  through  experimentation,  that  the  highest  order  G  that  this  program  can 
accommodate  is  Since  the  designed  FCS  G,jS  all  are  higher  order  than  8*^,  then 
order  reduction  must  be  accomplished.  The  method  used  to  reduced  the  order  of  the 
GijS  is  based  upon  obtaining  a  lower  order  approximation  from  the  Bode  plot  of  the 
designed  G',j  .  By  plotting  the  Lm(G'ij)s  on  a  Bode  plot,  a  straight  line  approximation 
can  be  made.  Using  this  approximation,  [-C,(jt)]m  =[G'„]mQtz(fc)  is  plotted  on  a 
Nichols  chart.  If  there  are  violations  of  the  3  dB  M-contour,  [G„]„i  is  modified  and  the 
modified  [Ci(k)]m  is  plotted.  This  method  is  repeated  until  a  satisfactory  lower  order 
\Gii]m  is  obtained.  Plots  of  for  all  FCS’s  are  shown  in  Figures  7.1  through 

7.8.  The  lower  order  GijS  are  shown  for  FCS^l  in  Equations  (7.1)  through  (7.3). 

Tho.se  for  FCS^2  are  shown  in  Equations  (7.4)  and  (7.5).  The  G,jS  for  FCS^^^21, 
FCS^22,  ans  FCS:?(^23  are  shown  in  Equations  (7.6),  (7.7),  (7.8),  respectively.  When 
these  reduced  order  G,^s  are  programmed  in  the  YF-16  simulation  program,  only  the 
pitch  channel  of  FCS#1  results  in  satisfactory  control  of  the  simulated  YF-16. 

(.3,00000-0.3)5''  -t- (3.18000-  0,3)5'* -<-(2.3.3940  -  03)52-f  (1.14730-03)5'  -|- (3.88580  -  04) 
^  57  +  (2450)5^  -I-  (4120000)55  (200000000)5'* 

(7.1) 
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r?i  — 


K'aalni  — 


[fj'salm 


( 1 .9200/}  +  07)5'^  +  (3,0720D  +  07)5^  +  (1.7280D  +  07)5^  +  (3.8400/)  +  06)S'  +  (2.4000D  +  05) 

5^  +  ( 1630)56  +  (688000)5-6  +  ( 19200000)5'' 

(7.2) 

(6.0000/)  -  01)5''  +  (6.6000/J-  01)5^  +  (8.9550/)  -  01)5-  +  (3.6390D  -  01)5'  +  (6  7035D  -  02) 
+  (2.7400D  +  03)5'’  +  (4.8280D  +  06)5^  +  ( 1.3888/)  +  09)5^  +  (4.8000D+  10)5'' 

(7.3) 

(7.56100  +  05)5''  +  (8.93760+06)5"  +  (1.49560  +  07)5"  +  (1.39500  +  07)5'  +  (9.11000+06) 
5"  +  (3.00000  +  02)5®  +  (6.00020  +  04)5®  +  (4.00020  +  06)5“' 

(7.4) 

(4.03860  +  0)5''  +  (5.89640  +  08)5"  +  (6.94240  +  08)5"  +  (3.95720  +  08)5'  +  (8.43750  +  07) 
5"  +  (2625)5®  +  (2790000)5®  +  (506250000)5“' 

(7.5) 

(1.80000+  10)5“'  +(5.40000+  10)5"  +  (5.85000  +  10)5"  +  (2.70000  +  10)5'  +(4.50000+09) 
5"  + (6030)5®  +  (9180200)5®  +  (271200000)5^  +  (1800000000)5" 
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(9.3750D+  15)5"  + (2.1562/)+  17)5"  +  (1.5000D  +  18)5'  +(2.8125D+  18)5' 

5"  +  (3.0 1.50D  + 04)5®  +  (3.0451/) +  08)5®  +  (1.0452/)+  12)5''  +  (1.5169D  +  14)5"  +  (5.6250D+  15)5" 

(7.8) 


-  400  -  3SO 


300  -250  -200  -  750  -  700  -SO 


[dcgj 

l^opjf7  70  .  PCS  if  7 


Figure  7.1.  FCS#1  [£i(o)]m 


7-2 


eo 


Phms^  fdcgj 
Loopjf2Sl  .  PCS  jf2 


Figure  7.2.  FCS#1  [£2(1)]^ 

1:3  Modeling  Problem 

With  the  [G.jim’s  implemented  in  the  simulation  program,  it  is  found  that  for 
the  pitch-up/pitch-down  maneuver,  q  output  has  an  unwanted  acceleration  compo¬ 
nent.  This  can  be  seen  in  Figure  7.9,  where  after  the  first  ramp-up/ramp-down  of 
the  response,  the  output  should  stay  at  zero  until  the  ramp-down/ramp-up  portion 
of  the  response  is  reached.  The  unwanted  acceleration  is  very  evident  in  Figure 
7.10,  where  there  is  a  zero  input  command.  Instead  of  staying  at  zero  there  is  an 
acceleration  that  causes  the  desired  output  to  be  displaced  downward.  The  cause 
of  this  unwanted  acceleration  is  traced  back  to  the  modeling  process.  As  stated  in 
Chapter  .3,  there  is  a  step  increase  in  the  thrust  at  the  start  of  the  simulation  that  is 
not  modeled.  During  system  identification  an  assumption  is  made  that  the  system 
under  consideration  is  at  steady  state.  This  assumption  is  made  to  insure  that  all 
components  of  the  output  are  due  only  to  the  declared  input.  In  this  case,  since  the 
step  increase  in  thrust  is  not  included  as  an  input  to  the  system,  the  plants  that  are 
identified  have  this  input  attributed  to  the  q  input.  The  thrust  input  is  attributed 
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Figure  7.3.  FCS#1  [>C3(2)]m 

to  the  q  input  because  it  results  in  an  output  about  the  y-body  axis  (q-axis).  Com¬ 
pounding  this  problem  is  the  mechanism  which  computes  the  elevator  trim  angle,  it 
does  not  account  for  the  added  control  moment  produced  by  the  thrust  vectoring. 
This  has  a  very  similar  effect  as  a  step  increase  in  thrust.  To  model  Vij(o)S  directly, 
the  u,j(o)S  are  required  to  be  equal  order  over  equal  order.  This  is  true  because 
MatriXx  will  not  accept  an  improper  transfer  function  and  the  u,j(o)S  are  to  be  in¬ 
verted  to  find  the  qij(o)S.  To  ensure  the  Utj(o)’s  are  equal  order  over  equal  order,  an 
extra  pole  is  inserted  at  —100  in  each  v,j(o).  The  assumption  is  that  this  pole  is  large 
enough  in  magnitude,  to  be  insignificant  in  the  u,j(o)S  and  would  be  removed  from 
the  qij{o)S.  As  it  turns  out  these  poles  are  significant  and  can  not  be  removed  from 
the  9ij(o)S  without  altering  the  frequency  response.  It  seems  that  the  model  requires 
at  lea.st  an  equal  order  over  equal  order  which  is  an  indication  that  there  are  other 
inputs  that  are  not  modeled.  According  to  Horowitz,  when  there  are  inputs  that  are 
not  modeled  one  would  expect  to  see  more  zeros  than  poles.  [9]  In  this  case,  the 
model  is  constrained  to  be  equal  order  over  equal  order  so  there  are  not  more  zeros 
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Figure  7.4.  FCS#2  [C^o)]m 

than  poles.  During  the  QFT  design  process,  these  plants  with  mismodeled  transfer 
functions  representing  the  pitch  channel,  cause  the  compensator  to  weight  the  con¬ 
tributions  of  the  pitch  input  incorrectly  and  when  the  compensator  is  installed  in 
the  simulation  program  the  unwanted  acceleration  is  produced. 
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Figure  7.5.  FCS#2  [£2(1)], 
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Figure  7.6.  FCS#21  [£3(2)]m 
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Figure  7.9.  Nonlinear  Simulation  Showing  Unwanted  Acceleration 
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7.11.  Nonlinear  Simulation  Showing  Unwanted  Acceleration  with  No  Input 
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Figure  7.12.  Nonlinear  Simulation  Showing  Corrected  Pitch- up/Pitch-down 
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ure  7.13.  Nonlinear  Simulation  Showing  Corrected  Pitch-up/Pitch-down 
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Figure  7.14.  Nonlinear  Simulation  Showing  Corrected  Pitch-up/Pitch-down 
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7.^  Conclusion 

This  chapter  describes  the  nonlinear  simulations  of  the  designed  PCS.  First  the 
reasons  for  and  the  methods  for  order  reduction  are  outlined.  Then,  the  prob¬ 
lem  with  modeling  encountered  during  the  nonlinear  simulation  of  the  manuevers  is 
described. 
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VIII.  Conclusion 


8.1  Discussion 

The  purpose  of  this  thesis  is  to  advance  the  state  of  the  art  of  FCS  design 
by  applying  nonlinear  QFT  to  the  design  of  a  FCS  to  control  maneuvers  of  an  ag¬ 
ile  aircraft.  To  do  this,  maneuvers  are  chosen  that  represent  maneuvers  that  are 
desirable  in  within  visual  range  combat.  There  are  two  maneuvers  specified,  a  pitch- 
up/pitch-down  maneuver  and  a  pitch-up/ velocity  vector  roll  maneuver.  By  using 
a  full  six-degrees-of-freedom  nonlinear  YF-16  simulation,  specifications  for  the  ma¬ 
neuvers  are  developed.  Once  the  specifications  are  developed,  a  set  of  maneuver 
simulations  are  flown  that  consists  of  enough  members  to  adequately  represent  the 
specifications.  Plant  transfer  functions  are  developed  through  the  use  of  system 
identification,  using  the  input/output  time  histories  of  the  maneuvers.  During  the 
system  identification  process,  the  qij(Q)  are  found  directly  and  are  constrained  to 
be  minimum-phase.  After  identifying  the  plants,  the  cascaded  MIMO  QFT  tech¬ 
nique  is  used  to  design  the  FCS.  Stability  boundaries  are  generated  through  the  use 
of  a  Matrix*  program.  Because  of  limited  algorithm  robustness,  it  is  found  that 
many  of  the  resulting  models  developed,  have  relatively  high  order.  When  the  effec¬ 
tive  plants  are  first  calculated,  they  contain  RHP  poles  that  make  it  impossible  to 
achieve  a  viable  design.  To  prevent  the  formation  of  these  RHP  poles,  F-boundaries 
are  developed.  When  the  loops  are  reshaped  using  the  F-boundries,  many  of  the 
extra  RHP  poles  disappear.  Before  shaping  the  loops  it  is  found  that  FCS  schedul¬ 
ing  is  needed  because  of  large  uncertainty  in  the  magnitude  and  phase  of  the  plants. 
During  loop  shaping  it  becomes  evident  that  the  design  will  not  met  the  specified 
because  of  multiple  large  RHP  poles  in  the  effective  plants.  Once  all  of  the  loops  are 
shaped,  prefilters  are  designed  using  a  modified  method  that  involves  plotting  all  of 
the  frequency  response  of  the  CLTF  and  comparing  them  with  the  desired  frequency 
response  model.  The  time  response  specifications  are  based  on  the  military  specifi- 
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cation,  ‘Flying  Qualities  of  Piloted  Aircraft’,  MIL-STD-1797A  [17].  Because  order 
of  the  designed  FCS  is  too  large  for  the  YF-16  simulation  program,  an  approximate 
FCS  is  developed  through  the  use  of  straight  line  approximations  on  the  Bode  plot. 
After  implementing  the  reduced  order  FCSs  in  the  nonlinear  simulation  program, 
problems  in  the  modeling  process  are  encountered  that  prevent  the  simulation  of  the 
FCSs. 

8.2  Conclusions 

Maxlike  The  system  identification  algorithm.  Maxlike,  in  Matrix*  is  very  useful  for 
the  problems  encountered  in  nonlinear  QFT.  Because  the  model  is  defined  by 
the  user,  maxlike  is  flexible  enough  to  give  transfer  function  plants  directly, 
that  are  required  for  the  nonlinear  QFT  design.  For  instance  the  inverse  of  the 
plant  matrix  is  identified  directly  in  this  thesis. 

Minimum-Phase  Plants  Using  Maxlike,  plants  that  appear  to  be  NMP  can  be 
modeled  as  MP.  The  validity  of  the  model  is  determined  by  the  application 
and  subsequent  simulations.  If  the  resulting  MP  transfer  functions  produce 
acceptable  fits  to  the  data  then  they  can  be  used  in  the  nonlinear  QFT  process. 

Stability  Boundaries  Computer  generated  stability  boundaries  are  very  useful 
during  the  QFT  design  process.  Using  computer  generated  stability  bound¬ 
aries,  the  process  of  plotting  and  cutting  out  templates  is  avoided.  By  letting 
Matrix*  generate  these  boundaries,  one  can  shape  the  loops  on  a  Nichols  chart 
with  the  boundaries  superimposed.  Care  must  be  taken,  when  using  these 
boundaries,  to  identify  the  actual  boundaries  of  plant  with  large  phase  uncer¬ 
tainty. 

F-Boundaries  Boundaries  can  be  developed  that  reduce  the  development  of  RHP  in 
effective  plants.  As  described  in  Chapter  4,  many  effective  plants  are  developed 
that  contain  large  RHP  poles.  Many  of  these  poles  are  prevented  by  redsigning 
the  controllers  using  F-boundaries. 
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MatriXx  There  is  a  limit  to  the  robustness  of  the  transfer  function  algorithms  in 
Matrix*  .  In  Chapter  4  an  experiment  is  described  that  shows  that  Matrix* 
algorithms  are  not  robust  enough  to  accommodate  the  manipulation  of  higher 
order  transfer  function. 

Model  All  inputs  to  the  system  under  consideration  must  be  accounted  for.  In 
this  case  a  step  input  of  thrust  was  not  properly  modeled  which  caused  an 
unwanted  ramp  in  the  output. 

8.3  Recommendations 

1.  System  identification  methods  tailored  to  QFT  MIMO  design  need  to  be  organ¬ 
ised  and/or  developed.  This  is  a  big  research  area  in  which  there  is  much  in  the 
literature.  There  are  many  CAD  packages  that  contain  system  identification 
algorithms,  two  available  at  AFIT  are  Matrix*  and  Pro-Matlab  [16]  [15].  Sys¬ 
tem  Identification  is  one  area  of  nonlinear  QFT  that  needs  to  be  developed.  If 
system  identification  methods  can  be  tailored  to  nonlinear  QFT  MIMO  design, 
then  the  time  required  for  a  design  can  be  shortened  considerably. 

2.  The  possibility,  that  plants,  that  are  previously  identified  as  NMP,  be  mod¬ 
eled  with  a  system  identification  method  that  can  constrain  the  model  to  be 
MP,  should  be  investigated.  In  some  cases  a  plant  modeled  as  NMP  may  be 
adequately  modeled  as  MP.  This  can  be  accomplished  by  using  system  identi¬ 
fication  with  the  constriant  that  the  identified  plants  be  MP.  If  the  results  of 
the  system  identification  justify  it,  the  plants  can  be  modeled  as  MP. 

3.  A  nonlinear  QFT  CAD  package  is  needed  that  can  overcome  limitations  that 
present  packages  have.  A  package  such  as  Matrix*  can  be  the  core  of  a  nonlinear 
QFT  CAD  package  if  robust  algorithms  can  be  developed  that  will  accommo¬ 
date  high  order  transfer  functions. 

4.  The  application  of  nonlinear  QFT  to  specific  maneuvers  seems  to  be  an  at¬ 
tractive  method  of  FCS  design  because  of  the  ability  to  rigorously  model  a 
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nonlinear  plant.  Another  attraction  of  this  method  is  the  method  of  collecting 
data  to  form  the  model.  The  data  can  be  collected  from  actual  flight  data 

5.  Thrust  vectoring  is  a  significant  control  in  operating  aircraft  at  high  angle 
of  attack.  Research  should  be  done  that  defines  the  boundaries  an  F-16  can 
operate  in  using  thrust  vectoring. 
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Appendix  A.  Example  Maneuvers 


This  appendix  contains  four  examples  of  the  input/output  time  histories  used 
to  develop  the  maneuvers. 
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Figure  A. 2.  Example  Maneuver  #111,  cont’d 
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Figure  A. 4.  Example  Maneuver  ^144 
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Figure  A. 5.  Example  Maneuver  #144,  cont’d 
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Figure  A. 6.  Example  Maneuver  ^^144,  cont’d 
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Figure  A. 8.  Example  Maneuver  #211,  cont’d 
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Figure  A. 11.  Example  Maneuver  ^243,  cont’d 
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Appendix  B.  System  Identification  Plots 


This  appendix  contains  plots  of  all  of  the  system  identifications  runs  where  ?/out 
is  the  output,  u  is  the  input,  and  u  is  the  input  based  on  the  identified  transfer 
function.  For  these  runs,  the  inverse  of  the  plant,  V,  is  the  form  of  the  transfer 
function  to  be  identified. 
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Figure  B.4.  u,  u,  and  pout  for  Plant  #1114 
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Appendix  C.  'Itansfer  Functions 


This  appendix  includes  Plants  for  all  Maneuvers  developed  through  system 
identification  in  Chapter  4. 
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-3.5242/? +04(s-  100)s^ 

(s  +  4.7048D  -  01)(s  +  2.2655D  -  01  ±  j4.8535L>  -  01)(s  -  4.3559D  -  01  ±  j3.1175D  -  01) 

1.1154£>+  12(s-  100)s‘’ 

(s  +  8.7164/?  -  01)(s  -  6.21 17D  -  02  d:  >1.0553/?  +  00)(s  -  8.9575/?  -  01  ±  >6.4361/?  -  01) 

-2.9600D  +  06(s  -  lOO)®** 

(s  -  1.2809D  -  01)(s  +  4.3613D  -  01  ±  >2.3457D  -  01)(s  -  9.3616D  -  01)(s  -  2.4845D  +  01) 

9.8882/?+  10(s-  100)s‘* 

(s  +  5.3386D  -  01)(s  -  4.6382/?  -  01  ±  >3.2645/?  -  01)(s  +  1.8024D  -  01  ±  >5.7787/?  -  01) 


<?120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 


C-2 


C.4  Plants  for  Maneuver  #2113 


9iio  - 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

a  5 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 


5.0933D+  10(s-  100)5“* 

(s  +  2.7955Z)  -  01)(s  +  9.9977 D  -  02  db  j2.7605D  -  01)(s  -  2.4760D  -  01  ±  jl.9164£)  -  01) 

2.0279^  +  07(5-100)5“* 

(s  +  2.3850D  -  01  ±  il.7510D  -  01)(5  -  8.2430D  -  01  ±  il.4965D  -  01)(s  +  4.6893D  +  00) 

2.4692D+  11(5-100)5^ 

(s  +  3.6370D  -  01)(s  -  3.2234D  -  01  ±  j2.3317D  -  01)(s  +  1.3121D  -  01  ±  j3.8847D  -  01) 

-2.9002Z}  +  04(s  -  100)5^ 

(s  +  1.3645D  -  01  ±  j3.9027Z?  -  01)(5  +  4.4281D  -  01)(s  -  3.8500L>  -  01  ±  j3.0683D  -  01) 

-2.8103D  -  01(s  -  6.5700D  -  02)(s  -  6.5700D  -  02)(s  -  6.5700D  -  02)(s  -  6.5700D  -  02)(s  -  100) 
(s  +  2.6435D  -  01  ±  j2.1 108D  -  01)(s  -  5.2433D  -  01)(s  +  1.4069£>  +  00)(s  -  2.1874I>  +  00) 

-7.3474D+04(5-  100)5“* 

(s  -  4.6061D  -  01  ±  i3.3534D  -  01)(s  +  2.0008D  -  01  ±  i5.4263D  -  01)(5  +  5.8344D  -  01) 

1.1691£)+  10(s-  100)5“* 

(s  +  2.9010D  -  01)(s  +  1,2295D  -  01  ±  j3.1325£)  -  01)(s  -  2.8235D  -  01  ±i2.1559£)  -  01) 

-1.5272£>  +  06(s-  100)5“* 

(s  -  2.4297D  -  01)(s  +  2.7085D  -  01  ±  j2.7A27D  -  01)(5  -  1.0468D  +  00)(s  -  9.3016D  +  00) 

3.6475D+  10(5-100)5^ 

(5  +  3.2153£>  -  01)(s  -  3.2165D  -  01  ±  j2.3318D  -  01)(s  +  1.3820D  -  01  ±  >4.0051^  -  01) 


Plants  for  Maneuver  #2114 

1.2987D+  12(s-  100)s“* 

(s  +  2.4479D  -  0l)(s  -  3.3362D  -  01  ±  j2.4190£>  -  01)(5  +  1.2852D  -  01  ±  jA.A2b2D  -  01) 

-5.4251£>  +  07(s-  100)s“* 

(s  +  1.4856D  -  01  ±  jl.4623Z?  -  01)(s  -  5.2524I>  -  01)(s  -  1.0916Z?  +  00)(s  -  9.9205D  +  00) 

1.2580D+  12(5-100)5“* 

(s  +  2.6912Z?  -  01)(s  -  3.4806D  -  01  ±  j2.5143D  -  01)(s  +  1.3591D  -  01  ±  j4.5319£>  -  01) 

-3.78720+04(5-  100)s“* 

(s  +  3.6930^  01  ±  j8. 13800  -  02)(s  -  3.89160  -  01)(5  -  1.230‘2O  -  01  ±  >3.74380  -  01) 

-3.45470  -  01(5  -  1.51140  -  04)(s-  1.5114O-04)(s-  1.51140  -  04)(5-  1.51140  -  04)(s  -  100) 
(s  +  2.6.5200  -  01)(s  +  1.02270  -  01  ±  >4.51580  -  01)(5  -  1.37180  +  00  ±  >2.44390  -  01) 

-3.86170  +  04(5-  100)s“* 

(s  +  3.86680  -  01  ±  >1.17340  -  01)(s  -  4.07810  -  01)(s  -  1.32160  -  01  ±  >3.93160  -  01) 


C-3 


?310  = 

9320  = 

9330  = 

C.6 

9iio  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.7 

9iio  = 

9120  = 

9130  = 


-7.5715D+  10(s-  100)s^ 

(s  +  4.7595/7  -  01)(s  -  1.7723D  -  01  db  j5.1192D  -  01)(s  -  5.4541D  -  01)(s  +  6.7049D  -  01) 

-9.5425D  +  05(s-  100)8^ 

(s  -  3.3190/7  -  01)(s  +  2.4448.0  -  01  ±  j2.4778D  -  01)(s  -  9.69390  -  01)(s  -  9.9078/7  +  00) 

-8.4457/7  +  10(s  -  100)s'‘ 

(s  +  5.3823/7  -  01)(s  -  1.8924/7  -  01  ±  j5.5047/7  -  01)(s  -  5.8244/7  -  01)(s  +  6.7985/7  -  01) 


Plants  for  Maneuver  ^1121 


1.49670+  ll(s-  100)s‘‘ 

Cs  -  1.34620  -  01  ±  j9.2905O  -  02)(s  +  1.69510  -  01)(s  +  4.82780  -  02  ±  il.66350  -  01) 

8.22550  +  07(8-  100)s'' 

(s  +  3.72980  -  Ol)(s  -  1.88450  -  01  ±  j3.92810  -  01)(s  -  4.61190  -  01)(8  -  3.36630  +  01) 

1.67680+  11(8-  100)8'* 

(s  -  1.20490  -  01  ±  j8.44890  -  02)(s  +  1.48040  -  01)(8  +  4.59290  -  02  ±  jl.45790  -  01) 

-1.93600+04(8-  100)8^ 

(s  -  1.03670  -  01  ±  j6. 36260  -  02)(s  +  1.85760  -  02  ±  jl.42370  -  01)(s  +  1.66370  -  01) 

-8.24400  -  01(s  -  5.86960  -  02)(s  -  5.86960  -  02)(s  -  5.86960  -  02)(s  -  5.86960  -  02)(s  -  100) 
(s  -  9.74260  -  02)(s  -  2.14090  -  01)(s  -  1.98320  -  02  ±  j5.26690  -  01)(s  +  2.66990+00) 

-1.53150+04(8-  100)8'* 

(s  -  9.59200  -  02  ±  j6. 15400  -  02)(s  +  2.27340  -  02  ±  jl.24910  -  01)(s  +  1.42530  -  01) 

1.86690+  ll(s-  100)s^ 

(s  -  2.57190  -  01)(s  -  1.02680  -  01  ±  j2.46910  -  01)(8  +  2.33850  -  01  ±  j2.0512O  -  01) 

3.75570+07(8-  100)s‘* 

(s  -  6.52340  -  01)(s  +  2.17810  -  02  ±  i7.0132O  -  01)(s  +  1.45900  +  00)(s  +  2.90670  +  01) 

4.39300+  10(s-  100)8“* 

(8  -  1.78600  -  01)(s  -  6.78260  -  02  ±  jl.71650  -  01)(s  +  1.56000  -  01  ±  jl.2606O  -  01) 

Plants  for  Maneuver  #2121 

1.04570+  10(s-  100)8“* 

^  -  3.80480-  01  ±  ;2.46710  -  01)(s  +  2.39100  -  02  ±  j4.8440O  -  01)(8  +  5.60820  -  01) 

9.65820+07(8-  100)8“* 

(8  +  4.12170  -  01)(s  +  7.93620  -  01  ±  j4.3930O  -  01)(s  -  3.73930  +  00)(8  +  3.76690  +  01) 

1.32780+  ll(s-  100)s“* 

(8  -  6.09900  -  01  ±  ;4. 19670  -  01)(8  +  8.64770  -  01)(s  +  3.73870  -  01  ±  i8.35270  -  01) 


C-4 


9210  = 

9220  = 

9230  = 

9310  = 

9220  = 

9330  = 

C.8 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 


4.8S55D  +  03(s-  100)s^ 

(s  -  3.0984£>  -  01  ±  j2.0883D  -  01)(s  +  5.8771D  -  02  ±  ;3.8353£)  -  01)(s  +  4.8641£>  -  01) 

-5.6243D-01(s-  lOO)®"* 

(s  +  2.8094D  -  01  ±  jl.9574D  -  01)(s  -  8  '■>167I>  -  Ui  ±  jQMMD  -  01)(s  +  4.7258D  +  00) 

-1.5323Z)+04(s-  100)s^ 

(s  -  5.4880D  -  01)(s  -  2.1775D  -  01  ±i5.2094D  -  01)(s  +  4.9626D  -  01  ±  j3.5676D  -  01) 

4.9337£)  +  09(s  -  100)s^ 

(s  -  2.55381)  -  01  ±  jl.6225D  -  01)(s  +  5.1558D  -  02  ±  j3.3548D  -  01)(s  +  4.0105D  -  01) 

-3.9581Ii  +  06(s-  100)s'* 

(s  +  5.0865D  -  01)(s  -  9.8382D  -  01)(s  +  1.G700D  +  00  ±  j2.3878D  -  01)(s  +  1.2210I>  +  01) 

-8.0893D  +  09(s  -  100)s^ 

(s  +  3.1565D  -  01  ±  jl.9473D  -  01)(s  -  3.7720D  -  01)(s  -  1 .3841Z)  -  01  ±  j3.6805D  -  01) 


Plants  for  Maneuver  ^^2122 

4.9833£)+09(s-  100)s’ 

(s  -  1.2290/?  -  01  ±  i8.5042D  -  02)(s  +  3.3480D  -  02  ±  jl.5659D  -  01)(s  +  1.7599D  -  01) 

-2.0187/?  +  06(s  -  100)s'‘ 

(s  +  5.3494D  -  02  ±  j 4.6898D  -  01)(s  +  5.7028/?  -  01)(s  +  5.0125/?  -  01  ±  j2.6525/?  +  00) 

_ -3.0669D-f  10(3-  100)s^ _ 

{s  -  2.7  137D  -  01)(s  -  9.3901/?  -  02  ±  i2.6477/?  -  01)(s  +  2.2985/?  -  01  ±  jl.6644D  -  01) 

1.1910D+04(s-  100)s‘* 

{s  -  4.5773/?  -  0l)(s  -  1.7627/?  -  01  ±  j4.3156D  -  01)(s  +  4.1912/?  -  01  ±  j2.7956D  -  01) 

-3.02590  -  01(s  -  3.5025D  -  02)(s  -  3.5025/?  -  02)(s  -  3.5025/?  -  n2)(s  -  3.5025D  -  02)(s  -  100) 
(.s  +  ?  4266D  -  01  ±  j2.0961D  -  01)(s  -  7.7948/?  -  01)(s  -  1.4887/?  +  00)(s  +  5.1650D  +  00) 

_ -5.0323D+04(s-  10P)s‘’ _ 

{s  +  7.8654D  -  01)(s  -  6.5125D  -  01  ±  >4.7236/?  -  01)(s  +  2.9557/?  -  01  ±  >7.7300/?  -  01) 

_ -1.6100D+  ll(s-  100)g‘’ _ 

{s  -  2.5813/?  -  01  ±  >5.7440D  -  01)(s  -  6.3431/?  -  01)(s  +  6.0260D  -  01  ±  >4.1533/?  -  01) 

1.1678/?+ 07(s-  100)s'' 

(.s  -  1 .5055D  -  01)(s  +  5.6659D  -  01  ±  >1.4565D  -  01)(s  -  1.6793/?  +  00)(f  +  2.9420/?  +  01) 

-3.8532D+  ll(s-  100)s^ 

(?  -  6.9427/?  -  01)(s  +  C.0581/?  -  01  ±  >3.7260/?  -  01)(s  -  2.5319D  -  01  ±  >6.6831D  -  01) 


C-5 


C.9  Plants  for  Maneuver  #2123 


9uo  = 

9120  = 

9130  = 

9210  = 

9220 

9230  = 

9310  = 

9320  = 

9330  = 

C.IO 

9iio  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 


1.0127D+  ll(s-  100)s'‘ 

(s  -  1.8705D  -  01)(s  -  7.7698D  -  02  ±  il.8049Z)  -  01)(s  +  1.6965D  -  01  ±  jl.4972£>  -  01) 

-3.9092£)+08(s-  100)s‘* 

(s  +  1.8317D  -  01  ±  jl.7932D  -  01)(s  -  7.1004D  -  01  ±  j3.8582D  -  01)(s  -  6.0134Z)  +  01) 

1.0475D+ll(s-100)s^ 

{s  -  1.8723D  -  01)(s  -  7.8462JD  -  02  ±  >1.80741)  -  01)(s  +  1.7049D  -  01  ±  >1.5106£)  -  01) 

-2.1273Z)+05(s-  100)s^ 

(s  +  6.8540D  -  01)(s  -  6.5060D  -  01  ±  >4.9194D  -  01)(s  +  3.4984Z)  -  01  ±  >8.0054£)  -  01) 

-2.4343D  -  01(s  -  3.4094D  -  02)(s  -  3.4094Z)  -  02)(s  -  3.4094D  -  02)(s  -  3.4094D  -  02)(s  -  100) 
(s  +  2.1599D  -  01  ±  >1.92000  -  01)(s  -  4.62220  -  01)(s  -  1.39680  +  00)(s  +  4.57960  +  00) 

-8.24630 +04(s-  100)s'‘ 

(s  +  4.94040  -  01)(s  -  4.94440  -  01  ±  >3.68490  -  01){s  +  2.70930  -  01  ±  >5.75670  -  01) 

5.00130+  10(s-  100)s^ 

(«  -  1.75910  -  01)(s  -  1.11410  -  01  ±  >1.82840  -  01)(s  +  1.90670  -  01  ±  >2.10250  -  01) 

-7.29740  +  06(5-  100)s^ 

(s  +  3.06790  -  01  ±  >2. 17260  -  01)(s  -  4.22580  -  01  ±  >8.43240  -  01)(s  -1.17120  +  01) 

5.01950+  10(s-  100)5“* 

(s  -  1.67890  -  01)(5  -  1.14310  -  01  ±  >1.78360  -  01)(5  +  1.89600  -  01  ±  >2.11030  -  01) 


Plants  for  Maneuver  #2124 

1.45800+  10(s-  100)5“* 

(5  -  1 .87300  -  02)(5  -  1 .05770  -  01  ±  >1.01590  -  01)(5  +  1 .13530  -  01  ±  >1.27830  -  01) 

-1.03900  +  08(5-  100)5“* 

(s  +  2.02740  -  01  ±  >1.87020  -  01)(5  -  7.09400  -  01)(5  -  1.43480  +  00)(s  -  2.20770  +  01) 

1.58620+  10(5-  100)5“* 

(5  -  3.74560  -  02)(5  -  1.00890  -  01  ±  >1.02420  -  01)(5  +  1.17910  -  01  ±  >1.30100  -  01) 

_ -1.73280+04(5-  100)5“* _ 

(s  +  1 .45670  -  01  ±  >4.42750  -  01)(5  -  3.98980  -  01  ±  >2.88610  -  01)(s  +  6.42830  -  01) 

-3.35870-01(5-  1.36440 -01)(5-  1.36440  -  01)(5  -  1.36440  -  01)(5  -  1.36440  -  01)(5  -  100) 
(s  +  2.51850  -  01  ±  >1.27450  -  01)(5  -  3.58600  -  01)(s  +  9.19400  -  01)(5  -  3.04590  +  00) 

-1.78130  +  04(5-  100)5“* 

(5+1 .4221 0  -  01  ±  >4.38350  -  01)(s  -  3.9181 0  -  01  ±  >2.82060  -  01)(s  +  6.45630  -  01) 


C-6 


9310  = 

9320  = 

9330  = 

C.ll 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.t2 

9110  = 

9120  = 

9130  = 


7.1611£>+09(s-  100)s‘‘ 

(s  -  3.2760D  -  02)(s  -  1.4603D  -  01  ± il.3766D  -  01)(s  +  1.5355D  -  01  ± il.9692I>  -  01) 

-4.2072D  +  06(s  -  100)s^ 

(s  +  3.0922D  -  01  ±  i2.9632D  -  01)(s  -  6.6991D  -  01)(s  -  3.0152D  +  00  ±  j2.1593D  +  00) 

7.1559£>+09(s-  100)s^ 

(s  -  3.1979D  -  02)(s  -  1.4481D  -  01  ±  jl.3512£>  -  01)(s  +  1.5208D  -  01  ±il.96eOZ:i  -  01) 


Plants  for  Maneuver 


6.9176D+  10(s-  100)s^ 

(s  -  1.3778Z)  -  01)(s  -  4.3416D  -  02  ±  jl.3077Z)  -  01)(s  +  1.0966D  -  01 .  9.1966/?  -  02) 

9.2104/?  +  07(s  -  lOO)®"* 

(s  +  1.0691D  -  01  ±  jl.8877/?  -  01)(s  -  3.1679/?  -  01  ±  j2.1725/?  -  01)(s  -  7.6254D  +  01) 

9.5546D+  10(s-  100)s^ 

(s  -  1.2543D  -  01)(s  -  4.5232D  -  02  ±  il.l805D  -  01)(s  +  1.0678/?  -  01  ±  j9.4032/?  -  02) 

-7.1332D  +  04(s-  100)s^ 

(s  +  1.2209/?  -  01)(s  -  3.2863/?  -  01)(s  -  6.8061/?  -  02  ±  j3.4784D  -  01)(s  +  3.7690/?  -  01) 

-8.5349/?  -  01(s  -  3.8906D  -  02)(s  -  3.8906/?  -  02)(s  -  3.8906/?  -  02)(s  -  3.8906/?  -  02)(s  -  100) 
(s  -  1.2964D  -  01)(s  -  4.2736/?  -  01)(s  +  8.1542/?  -  02  ±  j7.6122/?  -  01)(s  +  5.6592/?  +  00) 

-3.9110/? +  04(s-  100)s^ 

(s  +  3.3024/?  -  01  ±  J2.2082D  -  01)(s  -  4.1274D  -  01)(s  -  1.2529/?  -  01  ±  j3.9441/?  -  01) 

1.2138/?+  ll(s-  100)8“* 

(s  -  2.6774/?  -  01  ±  jl.8095/?  -  01)(s  +  3.2546/?  -  01)(s  +  9.6848/?  -  02  ±  ;3.5097/?  -  01) 

1.7441/? +08(s-  100)8“* 

(8  -  3.0282/?  -  01  ±  i5.7234D  -  01)(s  +  6.9995/?  -  01)(8  -  1.0075/?  +  01  ±  jl.l525/?  +  01) 

6.3922/?+  10(8-  100)8“* 

(8  +  2.4250/?  -  01)(s  -  2.0216/?  -  01  ±  jl.4028D  -  01)(8  +  7.6068/?  -  02  ±  ;2.4608/?  -  01) 


Plants  for  Maneuver  #2134 


8.6624D  +  09(8  -  100)8“* 

(s  -  1.5740D  -  01)(8  -  6.6605/?  -  02  ±  jl.51 18/?  -  01)(8  +  1.4330D  -  01  ±  jl.3204/?  -  01) 

-1.3972D+08(s-  100)s“* 

(s  +  2. 1103/?  -  01  ±  j  1.8624D  -  01)(8  -  7.8434D  -  01  ±  i4.8328/?  -  01)(8  -  5.0168/?  +  01) 

1.1021/?+  10(s-  100)8“* 

(.8  -  1.4740D  -  01)(s  -  6.3321/?  -  02  ±  jlAlOSD  -  01)(8  +  1.3537D  -  01  ±  il.2430/?  -  01) 


C-7 


9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.13 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 


9.8434D+03(s-  100)s^ 

(s  +  3.6872D  -  01  ±  j3.0014£>  -  01)(s  -  5.1884D  -  01)(s  -  1.5757D  -  01  ±  j5.0230D  -  01) 

-3.3909D  -  01(s  -  5.3926D  -  02)(s  -  5.3926£>  -  02)(s  -  5.3926D  -  02)(s  -  5.3926D  -  02)(s  -  100) 
(s  -  2.2842D  -  01)(s  +  1.8474D  -  01  ±  il.7786D  -  01)(s  -  1.9767D  +  00)(s  +  3.7740D  +  00) 

1.4643D+04(s-  100)s^ 

(s  +  3.9070D  -  01  ±  j3.1725£>  -  01)(s  -  5.3602D  -  01)(s  -  1.6671Z?  -  01  ±  j5.1710D  -  01) 

5.4599D  +  09(s-  100)s'‘ 

(s  -  1.0967D  -  01)(s  -  1.1005D  -  01  ±  jl.3887£)  -  01)(s  +  1.5777 D  -  01  ±  jl.8459£)  -  01) 

-4.1037IH- 06(s  -  100)s^ 

(s  +  2.4687D  -  01  ±  j2.2977D  -  01)(s  -  9.8941D  -  01  ±  j2.3693D  -  01)(s  -  7.2636D  +  00) 

8.4543Z)+09(s-  100)s^ 

(s  -  8.7034D  -  03)(s  -  1.4456Z)  -  01  ±  jl.2889D  -  01)(s  +  1.4180D  -  01  ±  il.8648D  -  01) 


Plants  for  Maneuver  ^^5^1141 

1.9415D+  ll(s-  100)s^ 

(s  -  2.4266D  -  01  ±  il.6285D  -  01)(s  +  2.9371D  -  01)(s  +  7.5519Z)  -  02  ±  j3.2063D  -  01) 

l,1435D  +  08(s-  100)5** 

(s  +  4.0762D  -  01)(s  -  2.4808Z)  -  01  ±  jA.9739D  -  01)(s  -  6.0617D  -  01)(s  -  1.3176D  +  02) 

3.8613D+  ll(s-  lOO)®** 

(s  +  2.6047D  -  01)(s  -  2.3558Z?  -  01  ±  jl.6108£»  -  01)(s  +  9.1054D  -  02  ±  j3.0329D  -  01) 

-1.1686£>+05(s-  100)s^ 

{s  -  3.6304D  -  03)(s  -  8.9804£)  -  01  ±  j7.3518£)  -  01)(s  +  OJOOIZ?  -  01  ±  jl.6560£)  +  00) 

-8.4623£>-01(s-  100)s‘' 

(«  +  4.0368£>  -  01)(s  -  3.5426D  -  01  ±  j5.9990£>  -  01)(s  -  7.6596D  -  01)(s  +  5.6342D  +  00) 

-1.2895£)+05(s-  100)5*’ 

(5  -  3.6414D  -  01)(5  +  3.4713Z?  -  01  ±  jl.4006D  -  01)(5  -  1.3580Z?  -  01  ±  j3.5770D  -  01) 

7.2517D+  10(5-  100)5*' 

(s  -  2.7972D  -  01  ±  jl.8797£>  -  01)(5  +  3.4702Z?  -  01)(5  +  9.3666Z?  -  02  ±  j3.6775D  -  01) 

2.7367D+07(s-  100)s'‘ 

(s  +  5.3956£)  -  01)(5  -  5.0012Z?  -  01  ±  j4.5010D  -  01)(5  -  4.4285D  +  00  ±  j5.6548D  +  00) 

7.9537i:i+  10(5-  100)5*’ 

(5  +  2.8597D  -  01)(s  -  2.4047D  -  01  ±  jl.6513D  -  01)(s  +  9.0345D  -  02  ±i2.9768£»  -  01) 


C-8 


C.14  Plants  for  Maneuver  #2144 


9110  — 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

a  15 

9110  = 

9120  = 

9l30  = 

9210  = 

9220  = 

9230  = 


1.7056D+  10(s-  100)s^ 

{s  -  3.6149D  -  01)(s  -  1.19190  -  01  ±  j3.4635D  -  01)(s  +  2.9122D  -  01  ±  j2.3130D  -  01) 

-1.3479D+07(s-  100)s^ 

(s  +  2.9166D  -  01  ±  jl.9896D  -  01)(s  -  3.1455D  -  01  ±  jl.0596D  +  00)(s  -  7.5903D  +  00) 

_ 4.2426D+  10(3-  100)3^* _ 

(3  -  3.7480Z?  -  01)(s  -  1.1982D  -  01  ±  j3.5607D  -  01)(s  +  2.9917D  -  01  ±  j2.4227D  -  01) 

8.5798D+ 03(3-  100)3^ 

(3  +  4.2718D  -  01  ±  j3.0871D  -  01)(3  -  5.7767D  -  01)(3  -  2.0096D  -  01  ±  j5.C228D  -  01) 

-3.8187D  -  01(3  -  3.6384Z)  -  02)(3  -  3.6384D  -  02)(3  -  3.6384D  -  02)(s  -  3.6384L>  -  02)(s  -  100) 
(s  -  1.7575D  -  01)(3  +  1.7196D  -  01  ±  il.7090£)  -  01)(3  -  1.7507Z)  +  00)(3  +  6.1784D  +  00) 

4.9964D  +  04(s-  100)3^ 

(s  +  4.9475D  -  01  ±  j3.9365D  -  01)(3  -  7.2197D  -  01)(s  -  2.4607D  -  01  ±  j6.9460D  -  01) 

1.0782D+  10(3-100)3“* 

(s  +  1.7219D  -  01)(3  -  2.0682D  -  01  ±  il.5389D  -  01)(3  +  1.1134D  -  01  ±  j2.4605Z)  -  01) 

-7.2518Z)  +  06(s-  100)3* 

(s  +  2.7400D  -  01  ±  jl.9853D  -  01)(s  -  2.4595D  -  01  ±  j6.3777D  -  01)(s  -  1.6295D  +  01) 

1.8414D+  10(3-100)3* 

(3  +  1.5304D  -  01)(s  -  1.9141D  -  01  ±  jl.4520Z)  -  01)(3  +  1.0656£»  -  01  ±  ;2.2413£>  -  01) 


Plants  for  Maneuver  #2211 


-2.95011) +  00(s-  100)s* 

(3  +  2.9358Z)  -  01)(3  +  1.5480£)  -  01  ±  j3.2042Z)  -  01)(3  -  1.2027Z)  +  00)(s  -  1.0018D+  01) 

5.4786D-01(s-  100)s* 

(s  +  5.3125D  -  01)(3  -  5.1484D  -  01  ±j2.S692D  -  01)(3  +  3.0772Z)  -  01  ±  j6.3867D  -  01) 

-2.7284£>-01(s-  100)s* 

is  +  3.1932D  -  01)(s  +  2.0343D  -  01  ±  j3.2116I>  -  01)(s  -  7.2340D  -  01  ±  >5.51541)  -  01) 

-1.6251D  +  00(3-  100)3* 

(s  +  1.3075D  -  01)(s  +  3.0966D  -  01)(s  -  4.9152D  -  01  ±  jl.7110D  -  01)(3  +  9.7227D+  00) 

-1.4213£)+00(s-  100)3* 

(s  -  2.4812Z)  -  02)(s  -  4.3214D  -  01  ±  >1.6076L>  -  01)(3  +  6.9982Z)  -  01  ±  >5.8600D  -  01) 

5.9419Z)-01(3-  100)3* 

(s  +  3.4379£>  -  01  ±  j2.8974D  -  01)(s  -  8.0989/)  -  01  ±  >1.4731D  -  01)(s  -  8.3791D  -  01) 


C-9 


9310  = 

9320  = 

9330  = 

C.16 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.I7 

9110  = 

9120  = 

9130  = 


-8.9112£>-01(s-  100)s'‘ 

(s  +  3.7669I>  -  01  ±  j2.2S68D  -  01)(s  -  1.0231D  +  00)(s  +  9.4073D  -  01  ±  jl.5585D  +  00) 

3.9907D-01(s-  100)s^ 

(s  -  6.13371)  -  01)(s  +  4.4891D  -  01  ±  j4.3733D  -  01)(s  -  3.0287£)  -  01  ±  i6.5823£>  -  01) 

-1.2938D-01(s-  100)s^ 

(s  +  3.4651D  -  01  ±  j2.2268D  -  01)(s  -  7.1375D  -  01)(s  +  6.4840D  -  02  ±  j7.1752D  -  01) 


Plants  for  Maneuver  #2212 


-2.5899D  +  00(s  -  100)s^ 

(s  +  1.5166D  -  01)(3  1  2.4715D  -  01  ±  i3.4495D  -  01)(s  -  1.3312D  +  00)(s  -  8A688D  +  00) 

9.89731>-01(s-  100)s‘‘ 

(s  +  4.1611D  -  01)(s  -  5.5316D  -  01  ±  j2.5872D  -  01)(s  +  2.7090D  -  01  ±  i8.1237D  -  01) 

-3.1298/) -01(s-  100)s^ 

(s  +  1.3388D  -  01  it  j4.5181D  -  01)(s  +  5.0277/)  -  01)(s  -  8.3874D  -  01  ±  i5.7037/)  -  01) 

-3.3921/) +  00(s-  100)s'‘ 

(s  -  2.6976/)  -  01)(s  +  2.7623D  -  01  ±  j3.1898/)  -  01)(s  -  8.3517/)  -  01)(s  +  1.9205/)  +  01) 

-9.2455/) -01(s-  100)8^ 

(s  +  3.9286/)  -  01)(s  -  6.0145D  -  01  ±  j3.4276/)  -  01)(s  +  2.6365/)  -  01  ±  >7.3212/)  -  01) 

1.1253/)  + 00(s-  100)5^ 

(8  -  3.3283/)  -  02  ±  i5.8464D  -  01)(8  +  6.2088D  -  01)(8  -  1.3179/)  +  00)(s  -  2.3133/)  +  00) 

-7.2723/) -01(8-  100)8^ 

(8  +  4.0604/)  -  01  ±  j3.9317D  -  01)(8  -  6.3725/)  -  01)(8  +  5.3457D  -  01  ±  >1.4046/)  +  00) 

4.7806D-01(8-  100)8^ 

(8  +  1.7766D  -  01)(8  +  1.3684/)  -  01  ±>5.3055/)  -  01)(8  -  5.6541/)  -  01  ±  >3.3676/)  -  01) 

-1.5160/) -01(8-  100)8“’ 

(8  -  3.6869/)  -  01  ±  >2.6408D  -  01)(s  +  5.5957/)  -  01)(8  +  4.9157/)  -  02  ±  >6.6485/)  -  01) 


Plants  for  Maneuver  #2213 


-2.4836D  +  00(s-  100)8“’ 

(8  +  1 .2980/)  -  01)(8  +  2.4176D  -  01  ±  >3.4535D  -  01)(8  -  1.4052D  +  00)(s  -  7.4094/)  +  00) 

1.8237/) +00(s-  100)8“’ 

(s  +  2.8421  D  -  01)(s  -  6.0401/)  -  01  ±  >1.4894/)  -  01)(8  +  5.7848/)  -  02  ±  >9.9416D  -  01) 

-3.3050D-01(8-  100)8“* 

(8  +  1.1351D  -  01  ±  >4.5473D  -  01)(8  +  5.0086/)  -  01)(8  -  8.3174/)  -  01  ±  >5.6753/)  -  01) 


C-10 


9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.18 

9no  = 

9120  = 

9l30  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 


-4.0889D  +  00{s-  100)s‘‘ 

(s  +  2.7891£>  -  01  ±  i3.3279D  -  01)(s  -  4.5009D  -  01)(s  -  8.0759£>  -  01)(s  +  2.0496Z?  +  01) 

-8.1639£>-01(s-  100)s'‘ 

{s  +  2.2160Z?  -  01)(s  -  6.1207Z)  -  01  ±  j2.SbS2D  -  01)(s  +  1.3478D  -  01  ±  j6.6846D  -  01) 

1.1498£>+00(s-  100)s^ 

(s  +  6.2583D  -  01)(s  -  8.8051Z)  -  02  ±  j6.2744D  -  01)(s  -  1.6361D  +  00  ±  j3.3446£)  -  01) 

-6.9867£>-01(s-  100)s^ 

(s  -  5.60961)  -  01)(s  +  5.0085D  -  01  ±  j5.2446D  -  01)(s  +  4.2774£>  -  01  ±  jl.3615D  +  00) 

7.2019D-01(s-  100)s'‘ 

(s  +  2.3748D  -  01)(s  +  1.9002D  -  02  ±  j6.2491i)  -  01)(s  -  8.3132D  -  01  ±  jl.8606D  -  01) 

-1.6282D-01(s-  100)s'‘ 

(s  +  6.6606D  -  01)(s  -  5.0197D  -  01  ±  i4.4413D  -  01)(s  +  9.6755D  -  02  ±  jT.lOOOZ)  -  01) 


Plants  for  Maneuver  #2221 

-1.9429£)  +  00(s-  100)s‘‘ 

(s  +  4.8464D  -  01)(s  +  6.5580D  -  02  ±  i4.8102D  -  01)(s  -  1.0495D  +  00)(s  -  6.8396D  +  00) 

1.1889£>+00(s-  100)s'* 

(s  -  8.7530D  -  02  ±  j5.1270D  -  01)(s  +  5.4148Z?  -  01  ±  >5.2092£>  -  01)(s  -  7.8964Z?  -  01) 

-4.9646Z)  -  01(5  -  1.9606D  -  04)(s  -  1.9606i)  -  04)(s  -  1.9606Z)  -  04)(s  -  1.9606Z)  -  04)(s  -  100) 

(s  -  5.2071Z)  -  01  ±  i3.9405D  -  01)(s  +  5.1564D  -  01  ±  i4.3131£>  -  01)(5  -  1.5883£>  +  00) 

-7.4804£>  +  00(s-  100)s^ 

(s  +  1.7820Z)  -  02  ±  j3.8243Z)  -  01)(s  +  4.1368D  -  01)(s  -  9.8501i)  -  01)(s  +  3.77750  +  01) 

-2.7257£)  +  01(s-  100)5“* 

(5  -  1.2744Z)  -  01  ±  j2.2370D  -  01)(5  -  1.2106Z)  +  00)(5  +  1.4476Z)  +  00)(s  +  4.1752D  +  00) 

1.0665D  +  00(s  -  2.3429D  -  04)(s  -  2.3429£>  -  04)(s  -  2.3429D  -  04)(s  -  2.3429D  -  04)(s  -  100) 
(5  +  4.6983Z)  -  01  ±  J3.9925D  -  01)(5  -  6.3426 £>  -  01  ±  j2.5682D  -  01)(5  -  2.6675D  +  00) 

-8.6710D  -  01(s  -  5.7000D  -  04)(5  -  5.7000D  -  04)(5  -  5.7000Z)  -  04)(s  -  5.7000Z)  -  04)(5  -  100) 

(5  -  4.5966D  -  01)(s  +  3.9878D  -  01  ±  j7.7992D  -  01){5  +  1.8366D  +  00  ±  j5.4867D  -  01) 

1.00047)  +  00(s  -  8.2390D  -  04)(s  -  8.2390D  -  04)(s  -  8.2390D  -  04)(s  -  8.23907)  -  04)(5  -  100) 

(s  +  7.56097)  -  01)(s  +  1.68707)  -  01  ±  j8.31 13D  -  01)(5  -  1.09687)  +  00  ±  j5.3846D  -  01) 

-2.83137)  -  01(5  -  6.00217)  -  04)(5  -  6.00217)  -  04)(5  -  6.00217)  -  04)(s  -  6.00217)  -  04)(5  -  100) 
(5  -  5.9992D  -  01  ±  >5.99957)  -  01)(s  +  9.84717)  -  01)(5  +  2.97517)  -  01  ±  >1.15017)  +  00) 


C-11 


C.19  Plants  for  Maneuver  ^2222 


9ilo  = 

9120  = 

9l30  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.20 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 


-1.8153D  +  00(s  -  100)s^ 

(s  +  5.4143D  -  01)(s  +  5.7060D  -  02  ±  j5.5880D  -  01)(s  -  9.5234D  -  01)(s  -  7.1895£>  +  00) 

3.0840D+00(s-  100)s^ 

(s  4-  6.7089D  -  03  ±  j4.7012D  -  01)(s  -  8.7928D  -  01)(s  +  6.6129D  -  01  ±  j8.1758D  -  01) 

-5.9350£)-01(s-  100)s^ 

(s  +  5.7275D  -  01  ±  i4.8794D  -  01)(s  -  5.5126D  -  01  ±  i5.1732D  -  01)(s  -  2.0234D  +  00) 

-1.1314D  4-  01(s  -  1.6440D  -  04)(s  -  1.6440D  -  04)(s  -  1.6440D  -  04)(s  -  1.6440D  -  04)(s  -  100) 
(s  4-  2.9263D  -  01)(s  4-  8.8330D  -  02  ±  j3.5541£>  -  01)(s  -  1.3611D  4-  00)(s  4-  4.8439D  4-  01) 

4.6105£)4-00(s-  100)s^ 

(s  4-  3.3577D  -  01)(s  -  6.9475D  -  01)(s  -  4.6729D  -  01  ±  j9.2394D  -  01)(s  4-  1.5515D  4-  00) 

9.7390Z?  -  01(s  -  8.3346D  -  04)(s  -  8.3346D  -  04)(s  -  8.3346D  -  04)(s  -  8.3346D  -  04)(s  -  100) 

(s  -  5.9821D  -  02)(s  4-  2.7710D  -  01  ±i2.1951£)  -  01)(s  -  1.1653D  4-  00)(s  -  2.1161D4-  00) 

-8.7081D  -  01(s  -  8.1004D  -  04)(s  -  8.1004D  -  04)(s  -  8.1004D  -  04)(s  -  8.1004D  -  04)(s  -  100) 
(s  -  4.5359D  -  01)(s  4-  3.7490D  -  01  ±  jSMlOD  -  01)(s  4-  1.6623D  4-  00)(s  4-  1.9444D  4-  00) 

2.4924jO  4-  00(s  -  3.0869D  -  04)(s  -  3.0869D  -  04)(s  -  3.0869D  -  04)(s  -  3.0869D  -  04)(s  -  100) 

(s  4-  5.6202D  -  01)(s  4-  2.4557D  -  01  ±  j9.0280D  -  01)(s  -  1.2320I>  4-  00  ±  ;4.6016£>  -  01) 

-3.0220D  -  01(s  -  8.5464Z?  -  04)(s  -  8.5464D  -  04)(s  -  8.5464D  -  04)(s  -  8.5464D  -  04)(s  -  100) 
(s  -  6.8546D  -  01  ±  i6.6713D  -  01)(s  4-  9.68t52Z)  -  01)(s  4-  2.1744D  -  01  ±  jfl.0850£)  4-  00) 


Plants  for  Maneuver  ^2223 

-1.7234D  4-  00(s  -  lOO)®"* 

(s  4-  5.2838D  -  01)(s  4-  5.3049D  -  02  ±  j5.5240D  -  01)(s  -  9.2115Z?  -  01)(s  -  7.0739D  4-  00) 

9.8879D  4- 00(s  -  100)5** 

(s  -  4.8765D  -  02  ±  i4.8385£>  -  01)(s  -  9.0839D  -  01)(s  4-  1.4761D  4-  00  ±  ;1.0698D  4-  00) 

-6.6338D  -  01(s  -  1.6962D  -  04)(s  -  1.6962Z?  -  04)(s  -  1.6962Z?  -  04)(5  -  1.6962D  -  04)(s  -  100) 
(s  4-  5.6377D  -  01  ±  j4.9117£)  -  01)(s  -  5.5530Z)  -  01  ±  j5,3013£)  -  01)(s  -  2.0673D  4-  00) 

-6.6391Z}4-00(s-  lOO)®** 

(s  -  1.8123D  -  01)(s  4-  2.8650D  -  01  ±  >3.5880D  -  01)(s  -  2.174411 4-  00)(s  4-  2.2337/1 4-  01) 

2A11GD  T-  00(5  -  7.SOG9iO  -  C4)(s  -  7.8559Z1  -  04)(s  -  7.8559Z?  -  04)(s  -  7.8559D  -  04)(s  -  100) 

(s  4-  2.9987£>  -  01)(s  -  3.5740D  -  01  ±  j7.0765£>  -  01)(s  -  8.4512D  -  01)(s  4-  1.9090D  4-  00) 

5.6123/1 -01(s-  100)5*’ 

(s  -  3.5992D  -  02  ±  i6.8443D  -  01)(5  4-  6.9204D  -  01)(5  -  1.0879/1 4-  00  ±  ;8.1422D  -  01) 


C-12 


93i0  = 

9320  = 

9330  = 

C.21 

9110  = 

9l20  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.22 

9110  = 

9120  = 

9130  = 


-8.4949D  -  01(s  -  1.9314D  -  03)(s  -  1.9314£>  -  03)(s  -  1.9314D  -  03)(s  -  1.9314D  -  03)(s  -  100) 
(s  -  4.3908D  -  01)(s  +  3.6764l>  -  01  ±  jQ.imSD  -  01)(s  +  1.5038D  +  00  ±  j4.8878D  -  01) 

3.8943D  +  01(s-  1.9607£>  -  04)(s  -  1.9607Z)  -  04)(s  -  1.9607D  -  04)(s  -  1.9607D  -  04)(s  -  100) 

(s  +  4.3615Z)  -  01)(s  -  6.8584D  -  01  ±  j6.5499D  -  01)(s  +  4.7205£>  +  00  ±  j8.m8D  -  01) 

-3.3028D  -  01(s  -  2.0247D  -  03)(s  -  2.0247D  -  03)(s  -  2.0247D  -  03)(s  -  2.0247D  -  03)(s  -  100) 
(s  +  9.3895D  -  01)(s  -  7.1954D  -  01  ±  >6.76837?  -  01)(s  +  1.7859D  -  01  ±  >1.0206D  +  00) 


Plants  for  Maneuver  ^22Z1 


-1.6659D  +  00(s-  100)s'‘ 

(s  +  5.5673D  -  01)(s  +  6.3690D  -  02  ±i5.6408£>  -  01)(s  -  2.1193D  +  00  ±  j8.1130£)  -  01) 

1.82777?  -  01(s  -  100)s^ 

(s  +  3.2424D  -  01  ±  j3.1969D  -  01)(s  -  1.46497)  -  01  ±  j5.4278D  -  01)(s  -  5.8115D  -  01) 

-2.0241D  +  00(s-  100)s‘‘ 

(s  -  5.53877?  -  01  ±  j5.0423D  -  01)(s  +  5.6086D  -  01  ±  >5.17037?  -  01)(s  -  9.1929D  +  00) 

4.1793D4-00(s-  100)s^ 

(s  +  9.3083D  -  02  ±  >4.35547?  -  01)(s  +  4.8704D  -  01)(s  -  2.75777?  +  00  ±  >3.29487?  +  00) 

-1.1283D-01(s-  100)5^ 

(s  -  8.4966D  -  02)(s  +  2.2962D  -  01)(s  +  1.48937?  -  02  ±  >3.09017?  -  01)(s  -  3.23717?  -  01) 

1.41817?  +  00(s  -  100)s^ 

{s  +  5.9598D  -  01  ±  >5.13360  -  01)(s  -  8.31260  -  01)(s  -  9.68490  -  01  ±  >1.09500  +  00) 

-1.70540  +  00(s  -  100)6“* 

(s  +  4.75910  -  01  ±  j3.6567D  -  01)(s  +  5.55030  -  03  ±  >6.86700  -  01)(s  +  1.06200  +  01) 

1.55360  -  01(6-  100)6“* 

(6  +  3.08140  -  01)(6  -  3.80320  -  01  ±  >1.21360  -  01)(6  +  4.57670  -  02  ±  >4.81030  -  01) 

-3.63500-01(6-  100)6^ 

(6  4  5.68460  -  01  ±  >4.84790  -  01)(6  -  7.99080  -  01)(6  -  6.24780  -  01  ±  >9.79500  -  01) 


Plants  for  Maneuver  #2232 


-1.76350  +  00(6-  100)6“* 

(6  +  5.44820  -  01)(6  +  8.33980  -  02  ±  >5.38690  -  01)(6  -  2.18190  +  00  ±  >1.34290  +  00) 

2.08000-01(6-  100)6“* 

(6  +  2.47330  -  01  ±  >3.14110  -  01)(6  -  1.31130  -  01  ±  >4.99230  -  01)(6  -  5.25760  -  01) 

-2.14820+00(6-  100)6“* 

(6  +  5.55860  -  01  ±  >5.21530  -  01)(6  -  5.94160  -  01  ±  >5.15350  -  01)(6  -  8.16450  +  00) 


C-13 


9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

933u  = 

C.23 

9110  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 


6.4792D+00(s-  100)s^ 

(s  +  8.5057£>  -  02  ±  j5.1043Z>  -  01)(s  +  5.5657D  -  01)(s  +  3.1834D  +  00  ±  i5.1336Z?  +  00) 

-8.1208D  -  02(s  -  100)s^ 

{s  +  2.5513D  -  01)(s  -  2.4239D  -  01  ±  jl.l628D  -  01)(s  +  1.5070D  -  02  ±  i3.2493I>  -  01) 

5.0114£>-01(s-  100)s^ 

(s  -  7.2056D  -  01)(s  +  6.2954D  -  01  ±  i5.3310D  -  01)(s  -  4.1589D  -  01  ±  jl.l3l3D  +  00) 

-1.9567D+00(s-  100)s^ 

(s  +  4.9735Z)  -  01  ±  j4.1525£)  -  01)(s  -  4.9046Z)  -  02  ±  j7.4676D  -  01)(s  +  1.0925Z?  +  01) 

2.3449D-01(s-  100)s^ 

(s  +  1.3332D  -  01  ±  jl.8243D  -  01)(s  -  1.7893D  -  01  ±  j5.0203D  -  01)(s  -  5.9489Z)  -  01) 

-3.7080D-01(s-  100)s'* 

(s  +  6.3141D  -  01  ±  i5.3781£>  -  01)(s  -  9.4434D  -  01)(s  -  6.0229Z)  -  01  ±  j9.9720D  -  01) 


Plants  for  Maneuver  #2233 


-2.0445D  +  00(s  -  100)s‘* 

(s  +  4.0216D  -  01)(s  +  7.4807D  -  02  ±i4.0424D  -  01)(s  -  1.3154P  +  00)(s  -  5.1161£>  +  00) 

5.50221>-01(s-  100)5^* 

(s  +  1.4265D  -  01)(s  -  3.2580D  -  01  ±  il.4363£>  -  01)(s  -  1.4795D  -  02  ±  jS.5365D  -  01) 

-2.3115Z?-l-00(s-  100)s‘’ 

{s  +  4.7614D  -  01  ±  j4,4189Z?  -  01)(s  -  5.2085£>  -  01  ±  j4.2885D  -  01)(s  -  4.6728£>  +  00) 

-1.6518Z?  +  01(s-  100)6“* 

(s  +  3.8575D  -  01)(s  +  7.1754D  -  02  ±  j4.3l90D  -  01)(s  +  2.6261D  +  00)(s  -  3.9335Z?  +  01) 

-6.5733£>-02(s-  100)6“* 

(6  +  2.9407D  -  0l)(s  +  9.9098D  -  02  ±  j3.9090D  -  01)(5  -  3.4419D  -  01  ±  j2.4331D  -  01) 

2.5835Z?-01(s-  100)s“* 

(6  -  4.6617D  -  01)(8  +  4.2696D  -  01  ±  i3.5597Z?  -  01)(6  -  3.9106£>  -  01  ±  j7.3085D  -  01) 

-2.0162I>+00(6-  100)8“* 

(8  +  3.i203D  -  01  ±  j2.978bD  -  01)(8  -  7.0477D  -  02  ±  j7.3070D  -  01)(8  +  6.9898Z1  +  00) 

6.09940-01(8-100)8“ 

(s  +  3.02640  -  01)(8  -  2.44020  -  02  db  j5.28530  -  01)(8  -  7.14050  -  01  ±  ;2.81990  -  01) 

-4.71340-01(8-  100)8“ 

(s  -  6.33070  -  01)(s  +  5.15760  -  01  ±  ;4.12460  -  01)(s  -  4.27140  -  01  ±  i7.46430  -  01) 


C-14 


C.24  Plants  for  Maneuver  #2241 


9iio  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330  = 

C.25 

9iio  = 

9120  = 

9l30  = 

9210  = 

9220  = 

9230  = 


-2.2838£>  +  00(s  -  100)s^ 

(s  +  5.4240D  -  01)(s  +  8.4502Z)  -  02  ±  j5.4508D  -  01)(s  -  1.2484£>  +  00)(s  -  5.5330D  +  00) 

4.1359D-01(s-  100)s‘‘ 

(s  +  2.5247D  -  01  ±  j3.5781D  -  01)(s  -  3.6976£)  -  01  ±  j6.7396D  -  01)(s  -  8.6179D  -  01) 

-1.1434D  +  01(s-  100)s‘’ 

(s  -  5.3552Z)  -  01  ±  J5.1478Z?  -  01)(s  +  5.5859D  -  01  ±  i5.2841D  -  01)(s  -  5.5250£l  +  01) 

1.5339D+01(s-  100)s^ 

(s  +  9.3075D  -  02  ±  j4.8177£)  -  01)(s  +  5.6289D  -  01)(s  -  2.3069D  +  00)(s  -  2.0711D  +  01) 

-1.8749D-01(s-  100)s‘' 

(s  +  2,3289£)  -  01  ±  j2A899D  -  01)(s  -  2.0736D  -  01  ±  i4.6028D  -  01)(s  -  6. HOOD  -  01) 

1.4223D+00(^  -  100)s‘' 

(s  +  7.1731Z)  -  01  ±  j5.8610£)  -  01)(s  -  6.9190D  -  01  ±  j9.1213D  -  01)(s  -  1.4384D  +  00) 

-3.3945£)  +  00(s  -  100)5“* 

(s  +  5.2628D  -  01  ±  i4.5915£>  -  01)(s  -  1.9237£)  -  01  ±  j7.9848£)  -  01)(s  +  1.4214D  +  01) 

3  3419£)-01(s-  100)5“* 

(s  +  2.8686 £)  -  01  ±  j3.4867D  -  01)(s  -  3.1628D  -  01  ±  i6.5375£)  -  01)(s  -  9.6510D  -  01) 

-4.1208Z)-01(s-  100)5“* 

(5  +  7.1871D  -  01  ±  i6.0533D  -  01)(s  -  6.2985D  -  01  ±  j9.7036D  -  01)(5  -  1.4006D  +  00) 

Plants  for  Manen  #2242 

-2.0808D  +  00(5-  100)5^ 

(5  4-  5.2832jD  -  01)(5  +  9.2051Z1  -  02  ±  j5.2770D  -  01)(s  -  1.1824D  +  00)(5  -  5.9138D  +  00) 

6.21131>-01(s-  100)5“* 

(s  +  1.8884D  -  01  ±  j3.3203D  -  01)(s  -  4.5737Z?  -  01  ±  j4.3776£»  -  01)(5  -  1.9648L>  +  00) 

1.85470  +  00(5-  100)s“* 

(5  -  5.3147D  -  01  ±  j4.92990  -  01)(5  +  5.47660  -  01  ±  j5.41270  -  01)(5  +  8.76120  +  00) 

3.06950+01(5-  100)5“* 

(s  +  7.39810  -  02  ±  i6. 17770  -  01)(s  +  6.88750  -  01)(5  -  3.40020  +  00  ±  jl. 06880  +  01) 

-1.91260-01(5-  100)5“* 

(s  +  2.35000  -  01  ±  >2.72320  -  01)(s  -  3.38750  -  01  ±  >4.68510  -  01)(5  -  9.74730  -  01) 

8.34050-01(5-  100)5“* 

(5  +  7.76830  -  01  ±  >6.50460  -  01)(5  -  7.19940  -  01  ±  >1.04390  +  00)(5  -  1.49200  +  00) 


C-15 


9310  = 

9320  = 

9330  = 

C.26 

9iio  = 

9120  = 

9130  = 

9210  = 

9220  = 

9230  = 

9310  = 

9320  = 

9330 


-3.0455D+00(s-  100)s^ 

(s  +  5.2659D  -  01  i:  j4.8198D  -  01)(s  -  2.8582D  -  01  ±  j8.7598D  -  01)(s  +  1.0278D  +  01) 

5.3681D-01(s-  100)s^ 

(s  +  2.5098D  -  01  ±  j3.4695D  -  01)(s  -  4.7938D  -  01  ±  j5.5161D  -  01)(s  -  1.7541D  +  00) 

-5.4511D-01(s-  100)s'‘ 

(s  +  7.3979Z)  -  01  ±  j6.3399D  -  01)(s  -  6.9379D  -  01  ±  j9.3224D  -  01)(s  -  1.7872I>  +  00) 


Plants  for  Maneuver  7^2243 


-2.0490D  +  00(s-  100)s^ 

(s  +  4.4250£>  -  01)(s  +  8.4514D  -  02  ±  j4.4206D  -  01)(s  -  1.0527D  +  00)(s  -  6.6044D  +  00) 

1.1573D+00(s-  100)s^ 

(s  +  8.0794D  -  02  ±  j2.9236D  -  01)(s  -  3.6507D  -  01  ±  i3.0765D  -  01)(s  -  2.8841D  +  00) 

1.3447Z?+00(s-  100)s^ 

(s  -  5, 1752D  -  01  ±  j4.5259D  -  01)(s  +  5.1816D  -  01  ±  jbA443D  -  01)(s  +  3.6340D  +  00) 

-1.1037D  +  01(s-  100)s'’ 

(s  +  2.0853D  -  02  ±  j5.5109D  -  01)(s  +  5.8328D  -  01)(s  +  3.4445/)  +  00)(s  -  1.4437/)  +  01) 

- 1.8376/) -01(s-  100)s'' 

(s  +  1.5931  D  -  01  ±  j2.4509D  -  01)(s  -  3.1062/)  -  01  ±  j2.9209D  -  01)(s  -  1.2132/)  +  00) 

3.6486D-01(s-  lOO)®'* 

(s  +  6.6048D  -  01  ±  j5.3166D  -  01)(s  -  9.5007D  -  01)(s  -  4.2695/)  -  01  ±  ;8.9941/)  -  01) 

-2.9339/)  +  00(s  -  100)s'’ 

(s  +  4.6804D  -  01  ±  ;4.3771D  -  01)(s  -  4.0220/)  -  01  ±  j8.3156D  -  01)(s  +  7.2709/)  +  00) 

1.0281/)+ 00(s-  lOO)®** 

(s  +  1.9814D  -  01  ±  j3.3401D  -  01)(s  -  4.3570/)  -  01  ±  j3.8085D  -  01)(s  -  2.9013/)  +  00) 

-6.5987D-01(s-  100)s‘’ 

(.s  +  7.4399D  -  01  ±  j6.1021/)  -  01)(s  -  5.5780/)  -  01  ±  ;9.5914/)  -  01)(s  -  1.3582/)  +  00) 


C-16 


Appendix  D.  Fall  Order  G’s 


This  appendix  includes  full  order  G’s  developed  in  Chapter  6  for  all  FCS’s 


G’s  for  FCS#1 


Sll(O) 


-{ 


1.2787D  -  03(s  -  5.1168D  -  01  ±  j  -  3.2726D  -  01)(j  -  2.7356£)  -  02  ±  j6.0847£)  -  01)  ^ 


_ (a  -  1.5000D  +  00)(a  -  5.0000D+  02) _ 

(a  -  1.0000£)  +  00)(j  -  1.0000D  +  02)(s-  1.8000D  +  02)(s  -  1.2000D  +  03  ±  jl.eOOOD  +  03) 

i(a-  1.2000D  +  03±i  -  1 .6000£)  +  03) ) 


S22( 


•>-( 


-1.3037D  +  10(s  +  4.7364D  -  01)(a  ~  6.8060D  -  01)(a  -  2.8766D  -  01  ±J  ~  7.6561D-  01)^ 


' _ (a  +  2.2818D  +  00)(»  -  5)(j  -  3.0821£>  +  02  ±  j4.3257£>  +  02) _ ) 

.  (a  -t-  8.1580D  -  0t)(s  -  7.Q000D  +  00)(s  -  l.OOOOD  +  Q2)(s  -  2.8207D  +  02  ±  j4.0795D  +  02) > 

(»  -  9.7937 D  +  02  ±  j  -  1.4039D  +  03) _ 


( 


(s  -  2.8207D  +  02  ±  j  -  4.0795D  +  02)(s  -  1.0071D  +  03  ±  j  -  1.4252D  +  03) 
( (s  -  4.0000D  +  03  ±  i9.1652D  +  03) ) 


\ 

/ 


333(2)  —  ^ 


( 


6.42160  -  01(s  -  5.42290  -  02  ±  j2.98250  -  02)(s  -  2.53590  -  01  ±  jl. 80710-  01)^ 


(s  -  3.00000  -  01  ±  j9.53940  -  01)(s  -  6.00000  -  01  ±  j  -  1.90790  +  00) 


(s  -  6.22730  -  02  ±j9.41810  -  02)(s  -  1.0000O+ 00)(5  -  2.25000  +  00  ±j  -  1.08970  +  00)/ 

{{a  -  1 .50000  +  01  )(a  -  9.00000  +  01)(a  -  6.2909O  +  02){a  -  8.57200  +  01  ±  jl .69930  +  03) \ 
I  (s  -  2.25000  +  00  + jl.0897O  + 00)(j- 5.50000  +  01)(3  -  1.00000  +  02)  ) 

' _ {a  -  8.57200  +  01  +  j  -  1.69930  +  03)(5  -  1.77800  +  03) _ ' 

,(a  -  3.08210  +  02  +  j  -  4.32570  +  02)(s  -  9.0001O  +  02)(s  -  1. 00720  +  03  +  jl. 42520  +  03) / 

1 


((*  -  1 


20000  +  03  +  j  -  1 .60000  +  03) 


+  03)) 


D-1 


G’s  for  FCS#2 


3ii(o) 


_  ^-7.5610£)  +  05(s-  3.00001)-  01  ±j9.5394D-01)(3  -  1.2027D  +  00)(s-  1.0018D  +  01)^ 


( (s  —  1.1 


1 


i-  1.7321D  +  02)) 


(s  -  1.0000D  +  02)(s-  1.0000D  +  02dt 


322(  1 ) 


2.6539D  +  08(s  -  2A966D  -  02)(s  -  4.3214D  -  01  ±  jl.6076D  -  01)(s  -  3.0000D  -  01  ±  ^9.53940  -  01) 


/(s-  3.0000D-  01  ±  j  -9.5394D-01)(a-  1.1859£)  +  00)(a-  1.5000D+  Ol)'^ 
\  (s  -  5.0000Z1-  01)(s-  1.2027£)  +  00)(s-  1.0003D  +  01)(s  -  1.0015D  +  01)  / 

"(3  -  1.2946D  +  01  ±j-  1.3229D  +  01)(»-  4.8346D+  02±  j  -  8.5657D  +  02)\ 
^  (s-3.5000D  +  01)(s-  2.2375D  +  02)(s-3.8752D  +  02±j8.2171£)  +  02)  / 


((s  -  1. 


(s-  1.2000D+03±j-9.0000D  +  02); 


G  for  FCS#21 


533(2) 


_  ^3.9150P  +  09(3  -  5.0197D  -  Ql  ±  j  -  4.4413D  -  01)(8  -  9.0000P  -  01  ±  j  -  4.35890  -  01)^ 


/  (3-4.50000  + 

\(s  -  1.0000O+  ( 

( 


00  ±  j2. 17940  +  00)(3  -  7.32390  +  00)(8  -  1.11680  +  02  ±  j  -  4.13430  +  01)  ^ 
00)(3  -  1.08330 +  01)(s-  5.63920  +  01)(3  -  8.33200  +  01  ±j  -  5.33140  + 01)  > 

_ (3-  1.11680  +  02  + j4. 13430  +  01)(3-  1.21540  +  02  +  jl.9180O  +  02) _ \ 

(s  -  8.3320O+  01  ±ji5.3314O  +  01)(s-  1.00000  +  02){s  -  7.85570  +  01  ±j  -  1.63100  +  02)/ 

1 


((3  -  1. 


(3-1 .20000  +  03  ±  jl .60000  +  03) 


+  03)) 


G  for  FCS#22 


/-i. 

.533(2)  =  I  - 


-1.09780  +  06(3  -  9.2200O  -  01  +  j  -  1.06180  -  01  )(3  -  5.42380  -  01  ±  jl  .07050  +  00) " 


/(3  -  4.50000  +  00  + >2.17950  +  00)(3  -  1.45980  +  01  >(3  -  2.20940  +  01  +  -  3.23870  +  00) 

\  (s  -  1.0000O  + 00)(s  -  1.19220  +  00)(3-  1.87580  +  01  ±j  -  4.35550  +  00)  / 

( _ (3  -  2.20940  +  01  ±  j3.23870  +  00)(3  -  3.41060  +  01) _ \ 

\(s  -  2.28480  +  01  )(s-  2.36440  +  01)(3-  1.0170O+  02)(3  -  3.00000  +  01  ±  jl.46970  +  02)/ 

((s  -  3.0000O  +  01  ±  J  -  1.46r7O  +  02)  ) 


D-2 


G  for  FCS#23 


a33(  2 ) 


/  -1.6589D  +  16(.?  -  2.2036D  -  02)(s  -  7.1634D  -  01)(i  -  3.0000D  -  01  ±  j  -  9.5394D  -  01) 

V 


/ _ (6-  -  7.4901P  +  00)(5-  3.0000D+  00±  j  -  9.5394P  +  00)(s  -  2.0000D+  01) _ 

V(s  -  2.5324D-  02)(s  -  1  .OOOOD  +  00)(4  -  1.0018D  +  01)(i  -  2.7825D  +  01  ±j  -  4.554l£)  +  01) 

/ _ (s  -  3.0000D  +  01  )(3  -  9.0740P  +  01  j:  jl  .61020  +  02) _ \ 

V(s  -  2.7825D  +  01  ±  J4.5541D+  01)(s-  7.1143D  +  01)(s-  6.3827D  +  01  ±  jl.6311D  +  02)/ 

((s  -  5.0000/9 +  03)(«-  1.0200D  +  04±  j  -  1.3600Z)  +  04) ) 


Appendix  E.  Prefilters 


This  appendix  includes  all  of  the  prefilters  used  in  this  thesis. 
F’s  for  FCS#1 


/ll(0) 


/22(1) 


/33(2)  — 


0.01(s  +  400) 
(s  +  4) 

0.1(s  +  50) 
(5  +  5) 

(0.05(5  +  80) 
(5  +  4) 


’s  for  FCS#2 


/u(o)  - 


0.05(5  +  100) 


/22(1)  — 


(5  +  5) 
5 

(5  +  5) 


F  for  FCS#21 


/33(2) 


(5  +  5) 


F  for  F(  'S#22 


/33(2) 


0.04(5  +  100) 

(«  +  i) 


F  for  I'(’S#2:i 
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